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In  this  paper,  we  study  certain  semidiscrete  methods  for  approximating  the 
solutions  of  initial  boundary  value  problems,  with  homogeneous  boundary  condi¬ 
tions,  for  certain  kinds  of  parabolic  equations.  These  semidiscrete  methods 
are  based  upon  the  availability  of  several  different  Galerkin-type  approxima¬ 
tion  methods  for  the  associated  elliptic  steady-state  problem.  The  properties 
required  of  the  spacial  discretization  methods  are  listed  and  estimates  of  the 
error  made  by  the  resulting  semidiscrete  approximations  and  of  its  time  deriva¬ 
tives  are  given.  In  particular,  estimates  are  given  that  require  only  weak 
smoothness  assumptions  on  the  initial  data.  Verifications  of  the  required 
properties  for  various  Galerkin-type  methods  are  also  given . 
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SIGNIFICANCE  AND  EXPLANATION 


Many  physical  situations  can  be  modelled  by  the  solutions  of  initial 
boundary  value  problems  for  parabolic  partial  differential  equations.  Examples 
of  such  situations  arise  in  the  theory  of  heat  conduction  and  other  diffusion 
processes.  The  physical  parameters  involved  are  often  dependent  on  the  time 
variable . 

The  construction  of  semidiscrete  approximations  to  the  solution  of  such 
parabolic  equations  is  studied  in  this  paper.  These  approximations  arise  from 
certain  spacial  discretization  techniques  and  they  are  governed  by  ordinary 
differential  equations.  Estimates  for  the  approximation  errors  are  given  for 
various  classes  of  initial  data,  including  classes  that  require  only  weak 
smoothness  assumptions. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


CONVERGENCE  ESTIMATES  FOR  SEMIDISCRETE 
PARABOLIC  EQUATION  APPROXIMATIONS 


Peter  H.  Sammon 

I .  Introduction 

Let  ft  be  a  bounded  domain  in  d-dimensional  Euclidean  space  with  a  sufficiently 

smooth  boundary  3ft  and  an  outward  pointing  normal  n(x)  =  (n,  ,  ■  •  •  ,  n .  Let 

—  1  d 

t  >  0.  We  shall  consider  semidiscrete  Galerkin-type  approximations  to  the  solution  of 
the  following  parabolic  initial  boundary  value  problem: 

d  d 

-u  =  L  (t)  u  =  -  y  D.  (a.  .(x,t)D.u)  +  ,)  a„ .  (x,t)D.u 

t  .  L.  ,  l  1]  3  Oi  1 

l , j=l  J  J  1  =  1 

(1.1) (i) 

+  aQ(x,t)u  in  ft  x  (0,t)  , 


under  one  of  the  following  boundary  conditions: 


(1.1) (ii) 


u(x,t) 


3ft 


0  or 


I 


ix  j=l 


a.  .  (x,  t) n. (x) D . u(x, t ) 
il  i  3 


3ft 


=  0, 


for  0  <  t  <  t 


and  with  the  following  initial  condition,  where  v  is  a  known  function: 

(1.1)  (iii)  u(x,0)  =  v(x)  for  x  e  ft 

(All  functions  considered  in  this  paper  will  be  real  valued) .  We  will  put  various 
kinds  of  restrictions  on  the  initial  data  function  v  later,  as  well  as  a  coercivity 


assumption  on  the  coefficients  of  L.  We  will  assume  that  {a..}  ,  { a„ . }  and 

13  .  .  ,  Oi  .  , 

1,3=1  1=1 

a„  are  sufficiently  smooth  functions  on  ft*  [0,xj,  that  a. .  =  a..  for  1  <  i,  j  <  d 

0  13  31  -  - 

d 


13  31 

form  a  uniformly  positive  definite  family  on 


and  that  the  matrices  [a. .] 

1]  i,3=l 

.  ■<  [ 0 , T 1 .  If  the  Neumann  boundary  conditions  are  under  consideration,  we  shall  also 
assume  that 

f  a..(x,t)  =  a(x,t)a..(x)  for  1  <  i,j  <  d  and  (x,t)  e  ii  x  [0,r]  , 

13  13  - 


(1.2) 


|  aoi  (x'fc) 


0  for  0  <  t  <  t  , 


3ft 
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where  [a^..  (x) 


is  a  symmetric  matrix  of  sufficiently  smooth  functions  on  f! 


i.  j=l 

and  a(x,t)  is  a  sufficiently  smooth  function  on  B  x  [ 0 , T ] . 

Suppose  we  consider  the  following  elliptic  boundary  value  problem  associated 
with  (1.1):  Given  0  t  <_  t  and  a  suitable  function  f,  find  a  function  w  which 
satisfies  L(t)w  =  Lw  =  f  on  fl  and  the  appropriate  homogeneous  boundary  condi¬ 
tions  on  3Q.  There  are  many  well  studied  techniques  for  finding  an  approximation, 
in  a  finite  element  space,  to  the  solution  of  this  problem  (see  the  surveys  in  [2] 
or  (3),  for  instance).  Moreover,  the  authors  of  (3)  have  shown  that  it  is  possible 
to  take  such  a  technique  and  use  it  to  generate  a  time  continuous  family  of  approxima¬ 
tions  to  a  solution  of  (1.1),  at  least  if  L  has  time  independent  coefficients. 

They  proved  that  this  semidiscrete  approximation  to  (1.1)  can  be  a  good  one  at 
positive  times,  even  if  the  initial  data  function  v  is  not  smooth  on  9.  or  does 
not  satisfy  suitable  boundary  compatibility  conditions.  This  indicates  that  the 
smoothing  property  of  the  parabolic  problem  can  be  utilized.  (It  is  known  that 
solutions  of  (1.1)  are  smooth  in  space  and  time  for  positive  times,  even  if  v  is 
not  smooth  and  compatible  in  space).  They  also  proved  uniform  in  time  results  for 
naturally  restricted  classes  of  initial  data. 

We  intend  to  study  the  case  of  time  dependent  non-selfadjoint  operators  in  this 
paper  and  show  that  many  of  the  results  of  (3)  are  still  valid,  under  similar 
hypotheses  on  the  data  of  the  problem.  We  will  begin  by  setting  some  notation  and 
then  we  will  describe  some  results  concerning  the  smoothing  properties  of  the  problem 
given  by  (1.1).  (These  results  are  known,  but  we  shall  give  some  derivations  that 
will  allow  us  to  prove  similar  results  for  the  discrete  setting).  We  will  then  make 
some  abstract  hypotheses  concerning  finite  element  approximations  to  the  associated 
elliptic  problem,  define  the  semidiscrote  approximation  and  prove  (optimal  order) 


convergence  results  for  restricted  classes  of  smooth  initial  data  as  well  as  for 
nonsmooth  initial  data.  It  will  then  be  shown  that  many  of  the  known  Galerkin-type 


approximations  satisfy  the  abstract  conditions.  All  of  these  estimates  will  be  done 
in  a  l?  (Q) -setting  but  we  will  conclude  by  discussing  some  error  estimates  in  the 
maximum  norm. 

We  refer  the  reader  to  [3]  for  a  discussion  of  related  work  done  by  other 
authors.  This  work  represents  an  extension  of  work  done  in  [8]  under  the  supervision 
of  Professor  J.  H.  Bramble. 

We  conclude  this  Section  with  an  observation  concerning  scaling  arguments.  If 

-Kt 

u(t)  is  a  solution  of  (1.1)  then  w(t)  e  u(t)  is  a  solution  of  the  following 
evolution  equation: 

w  (t)  +  { L ( t )  +  K)w(t)  =  0  , 

for  any  K  <=  TR.  This  relation,  as  well  as  a  similar  one  which  will  hold  for  the 
semidiscrete  approximation,  will  be  used  later. 

We  will  use  the  symbol  C  to  denote  a  generic  positive  constant  throughout  this 

work  and  we  will  define  /  (•)  '  0  if  m  <  m,  . 

.  2  1 


II .  Parabolic  Regularity 

We  will  now  fix  some  notation  and  discuss  some  of  the  properties  of  the  paraboli 
problem  defined  by  (1.1).  We  will  not  cite  specific  references  for  the  results  con¬ 
cerning  the  elliptic  equation  theory,  but  most  of  the  statements  we  make  may  be  found 
for  instance,  in  [6]. 

We  let  W^'15  5  W^'^(P.)  be  the  usual  L^(f2) -based  Sobolev  spaces  on  whore 

1  p  <  «°  and  l  >_  0  is  integral.  We  give  them  their  usual  norms,  denoted  by 

H  Si  ,2  2 

11*11  .  We  will  write  H  for  W  (the  L  (Q)-based  analysis  will  play  the 

»  P 

largest  role  in  our  work),  II  *  II  ^  for  II  *11  ^  ^  and  II  'II  for  II  *11^.  We  will  also  let 

2  2  1 
11*11  denote  the  L  (P.)  -*•  L  (0)  operator  norm.  We  will  use  the  space  H^,  which  is 

the  subspace  of  H'*'  consisting  of  functions  that  vanish  (in  the  sense  of  trace)  on 

3Q,  if  we  are  considering  the  Dirichlet  boundary  conditions.  We  will  also  need  to 

2  1 
consider  the  space  L  (30),  whose  norm  will  be  denoted  by  II  *11  and  H  (3?.), 

0 1  o 

whose  norm  will  be  denoted  by  11*11  .  We  let  (*,*)  and  <*,*>  denote  the  usual 

2  2 

L  (0)  and  L  (30)  inner  products,  respectively.  Finally,  if  X  is  a  Banach  space, 

£*f«r  £  +  e 

we  will  let  C  (|a,b],X)  and  C  ((a,b],X)  denote  the  usual  spaces  of  X-valued 

functions  that  have  a  Hblder  continuous  £-th  derivative,  where  t  >  0  is  integral 

and  0  ^  c  <  1  and  B(X)  denote  the  usual  Banach  space  of  bounded  operators  on  X. 

2 

We  shall  assume  throughout  this  work  that  the  initial  data  function  v  £  L  (0). 
We  now  need  some  concepts  that  pertain  to  elliptic  equation  problems  that  are 
associated  with  (1.1).  We  begin  by  defining  some  bilinear  forms  on  *  H^.  Let 
0  <_  t  <_  t  and  set 

d  d 

D  (t  J  (*  ,  * )  5  7  (a.  ,D.  (*)  ,D.  <• ))  +  7  (a.  D.  (•),(•))  +  (a  (•),(•)), 

...  in  1  i  01  i  0 

i,]=l  J  J  1=1 

d  d 

D*  ( t)  ( •  ,  • )  I  y  (a,  ,D.(*),D.  (•))  -  y  (a  ,  D.  (■),(•) ) 
i,j=l  13  ^  i  i=l  01  1 

+  ((ao-  .f  Dia0i)(-,-(*))  ■ 

1=1 


-4- 


-iw„  ,-w,  ’WW.  w*' 


As  usual,  D(t)  is 

for  L(t). 

D*  (t) 

bears 

defined  by 

c 

* 

ft 

=  L  ( t) 

L (t)  (•)  5 


associated  with  a  weak  formulation  of  a  boundary  value  urob 1  <.■ 
a  similar  relation  to  L*(t),  the  formal  adjoint  of  Lit), 

-  G(t)  where 


d  d 

G(t)(-)  H  I  aoiD.(-)  +  -  (  l  D  a  )(•)  . 

i=l  i=l 

1  vd 

We  shall  assume  that  a^  -  (a^  -  —  £  D^a  )  —  ^  °n  ^  x  10>T1  if  we  are  working 

i=l 

with  the  Dirichlet  boundary  conditions  or  that  aQ  >  0  on  2  *  [0,t]  if  we  are  vrork- 
int  with  the  Neumann  boundary  conditions.  Thus  D(t)  and  D*(t)  are  (strongly) 
coercive  forms  over  if  we  have  the  Dirichlet  boundary  conditions  or  over  H ^  if 

we  have  the  Neumann  boundary  conditions. 

2 

Since  we  wish  to  regard  L,  L*  and  L  as  unbounded  operators  on  L  (.1)  ,  we 

must  discuss  their  (common)  domain  of  definition  D  .  If  the  Dirichlet  conditions  are 

L 

2  1 

under  consideration,  we  let  D  =  H  n  H  .  If  we  are  considering  the  Neumann  problem, 

L  o 


2  2  _ 
we  let  D_  =  H  n  {w  e  H  :  )  a. .n.D.w!  =  0}.  Then  L,  L*  and  L  are  indeed 

L  i,j=i  13  1  3  Ian 

2  2 
closed,  unbounded  operators  on  L  (f!)  with  a  common  domain  D^,  L*  is  the  L  CD- 

adjoint  of  L  and  L  is  selfadjoint.  We  will  give  the  II  *11  ^-norm  and,  for 

convenience,  we  will  give  G(t)  and  G*(t)  5  -G(t)  the  domain  D^. 

2 

We  will  now  identify  a  space  intermediate  to  L  (SI)  and  that  will  be  useful 

later.  Let  0  <  t  <  x.  Since  L(t)  is  selfadjoint  and  positive  definite,  we  can 

-  i  -  1  2 

use  spectral  theory  to  define  L  (t)  on  its  domain  D(L‘(t))  e  L  (f!)  and  we  can  give 

_  i 

the  latter  the  norm  II  L  ‘  (t  )  (•  )ll  .  If  g  e  D^,  then 

II  L  "  gll  2  =  (Lg,g)  =  j  (D  +  D*)  (g,g)  =  II  gll  2 


where  we  have  used  -  to  denote  a  norm  equivalence.  Thus  it  can  be  seen  that 


-5- 


D(L'it))  is  the  II  *11  ^-norm  closure  of  and  the  II  L  ' ( t )  ( •  )ll  -norm  is  equivalent  to 

the  II  *11  -norm.  We  will  write  to  denote  this  space,  which  is  H1  if  the 

1  *  o 

Dirichlet  boundary  conditions  are  under  consideration  or  if  the  Neumann  conditions 

are  being  used  and  we  will  give  it  the  II  *11  ^-norm. 

We  let  T ( t )  denote  the  solution  operator  for  the  elliptic  boundary  value  problem 

associated  with  L(t),  for  0  <  t  <  T.  Thus  T(t)  is  an  operatt  r  from  right 

i.  1+2 

hand  sides  g  e  H  ,  for  any  !  0,  to  solutions  in  H  n  that  satisfy 

L(t)(T(t)g)  =  g.  We  define  T*(t)  and  Tit!  analogously  and  note  that  T* (t)  is 
2 

indeed  the  L  (fi)-adjoint  of  Tit)  and  that  Tit)  is  selfadjoint.  (We  will  continue 
to  use  the  symbol  *  to  denote  adjoints  taken  with  respect  to  the  L2(Sl)-inner  product). 

Let  j  0  and  let  L^'(t)  =  (— )  Lit)  denote  the  operator  obtained  by  dif¬ 
ferentiating  the  coefficients  of  Lit)  with  respect  to  time.  We  give  this  family  of 
operators  the  domain  Dl  and  we  define  L*^*(t)  and  L^*(t)  similarly.  If  we 

regard  Tit)  :  H'  ■+  H?'+2  n  Dr  for  some  i  >  0,  we  can  verify  that  T(1)  -  (~)T 

1  —  dt 

exists  (in  the  operator  norm)  and  is,  in  fact,  T(1)  =  -TL(1)T.  We  can  continue 
differentiating  and  show  that  T<:il  it)  exists  for  each  j  >_  0.  Similar  statements 
hold  for  T*(t)  and  Tit). 

We  can  now  use  the  work  of  Sobolevskii  [10)  to  study  the  solution  of  (1.1).  The 

solution  u  of  (1.1)  can  be  described  via  a  family  of  fundamental  solution 
2 

operators  U(t,s)  r  B  (L  P)  )  ,  defined  for  0  <_  s  <_  t  <  t  .  In  fact,  the  operators  U(t,s) 

2 

are  stronglv  continuous  in  L  P)  for  0  <  s  <  t  <  are  continuously  differentiable 

2 

in  each  variable  in  B  ( L  if! ) )  and  have  range  in  for  0  <_  s  <  t  <  T  and  are 

characterized  by  the  following  equations: 

(2.1)  U(t,s)  +  L(t)U(t,s)  =  0  and  U(s,s)  =  I  for  0  <_  s  <  t  r 

The  unique  solution  of  (1.1)  is  given  by  uit)  =  U(t,0)v,  for  0  ^  t  <_  x.  We  note 
that  U(t,;)U(:,s)  =  U  ( t ,  s)  for  1  r  s  c  <t<T  and  that  —  U(t,s)  =  U(t,s)L(s) 

—  —  d  S 

on  Dt  for  0  ^  s  "  t  ^  . 


-G- 


We  can  use  such  fundamental  solution  operators  to  describe  the  solution  of  a  more 
general  problem  than  (1.1).  The  following  equations: 

(2.2)  (i)  wt  +  L(t)w  =  f(t)  for  0  <  t  <_  t  and  w(0)  =  wQ  , 

where  wQ  e  L2  (D )  and  f  e  C*;  (  [0,t  1 ,1?  (Q) )  are  known  and  >  0,  have  a  unique 
solution 

(2.2)  (ii)  w  t  C(  [0,xl  ,L2C!)  )  n  C1  ( (0 ,  t  1  ,1?  (fl> )  nC((0,T),DL> 


given  by 


t 


(2.  3) 

w(t)  = 

U(t,0)w  +  j  U(t,s)f(s)ds  for  0  <_  t 

0 

i  T 

If  wQ  • 

V 

then  w  ■_ 

C([0,t),Dl)  and  if  w^  =  f ( 0 )  =  0,  w 

r  C  1  ( ( 0 ,  t  ]  ,  D  )  for 
L 

some  e, 

■'  0. 

Moreover , 

)  +  1  2 

if  ?.  >  0  and  f  c  C  ( 10,  t  )  ,L  (.1)  )  , 

then 

w 

( [0,T] 

l,L2C.)) 

and 

(2.4) 

;.+i 

|  w  ( t  >  II 

<■  CO)  t  '  1  (II wJ  +  sun  II  (-p-)  f(s)ll) 

—  0  ds' 

O'  s<t 

for  o  <•  t  T 

for  some  constant  C(9.). 


(We  will  often  write 


w(j) (t) 


ut;i>(t)  (— 1  3 u ( t)  and  U(^(t,s)  f--)  U(t,s)  for  j  >  0,  in  the  future). 

'■dt'  jv  " 

We  also  note  that  other  results  from  (10)  show  that  if  w  has  the  continuity 

properties  described  in  (2.2)(ii)  and  w^_  +  Lw  =  0  for  0  <  t  i  where 

w  (0 )  -  hJ  =  D  (I.  ‘(0) ) ,  then  w  <  c'  (  [0,  t  ]  ,  L2  (3)  )  and  (tw(t))  *:  C*  ( [0, t  3  ,  D^)  ,  for 

some  <■  >  0.  These  results  allow  us  to  derive  the  following: 

2 

Proposition  (2.1):  If  w^  •  D^,  f  ♦  C  (10#t),L  (.•!))  for  some  •  >  0  and 

L(0)w  -  f  <  0)  r  ,  then  the  solution  w  of  (2.2)  is  in  C  ^((0,t],D  )  for  some 

0  * 


Proof:  Let  •/  +  L(t)y  =  0  for  o  -  t 

z  -  w(t)  -  w()  +  ty(t)  for  0  ^  t  _  :  , 


r  and  y(0)  =  LfOfw^ 


f  (0) 


H^.  Then  if 


,  i  v:  -v. 1 1 


there  is  an 


2 


0  so  that 


zr  +  Lz  =  (f  -  Lw  +  y  +  (L  -  L)  (ty) )  5  g  e  C  2  (  [0,T  )  ,l/  <  .)  j 
*  3 

urvi  z(0)  =  =  0.  Thus  z  •  C  (I0,t),D  )  for  some  •  :>  0  and  the  re  -jit 

L  3 

Wo  will  soon  want  to  examine  further  conditions  under  which  the  solution  of 
is  smooth  at  t  =  0  and  we  will  al  -  want  to  estimate  the  size  of  this  sol  itior.  :s 
terms  of  data.  We  shall  study  another  result  about  solutions  of  (2.2)  to  carry  out 
this  analysis. 

1  2 

Proposition  (2.2):  If  f  ■  C  ([0,t],L  (11))  is  such  that  there  is  a  (unique)  soiutio 
1  2 

w  e  C  (|0,tJ,L  (T. ) )  n  C([0,t],D  )  of  (2.2)  satisfying  w(0)  =  w„=  0,  then  for  each 

n  0 

•  >  0,  there  is  a  constant  C  =  C(e)  >  0  so  that 


(2.5) 

II  w  ( t )  II 

<_  C  sup 

II  Tf  (s)ll  +  e  sup  s  II (Tf )  <s)ll 

for 

0<_s<t 

0<s<t  " 

(2.6) 

llw  (till 

C  sup 
0<  s<  t 

llTf(s)ll  +  «  sup  sll  Tf  (sill 

1  0<s<  t  s  1 

for 

P roo f :  We  see  that  Tw^  +  w  =  Tf  and  w(0)  -  0.  We  will  analyze  this  equation  usin 
energy  techniques  that  were  used  in  a  similar  argument  in  [3). 

We  first  see  that  if  0  <  t  <  i,  then 


(2.7) 


,  t  t 

i2_f  mii2.  r  _  d 


til  w ( t ) II  =  /  II  wfl  ds  +  /  s  —  llw(s)l|  ds 
'  J  ds 


The  equation,  (2.2),  shows  that 

(2.8) 


7  -r-  II  vHI  2  <  q-  -p-  II  wll  2  +  (Tw  ,w  )  =  (Tf  ,w  ) 
2  ds  —  2  ds  s  s  s 


=  (Tf ,w)  -  ( (Tf )  ,w)  for  0  <  s  < 

s  s  — 


Thus,  if  we  integrate  (2.8),  we  find  that 


(2.9) 


/  s  ^-llw||2ds  <_  jll  wll  2  +  CtllTfll2  +  C  /  II  Tfll  2ds  +  C  /  II  wll  2ds 


+  q-  t  sup  s2  II (Tf )  II 
2  0<s<t  s 


Since  Tw^  +  w  =  Tf  -  TGw,  we  see  that 
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(2.10)  (Tw  ,w)  +  llwll  2  =  i  (Tw,w)  -  ^  (T!1)w,w)  +  II w!l  2 

s  2  s  2 

=  (Tf,w)  -  (TGw,  w)  for  0  <  s  <_  t 

We  now  analyze  the  individual  terms  in  (2.10).  We  have  the  following  estimates,  where 


(T^’w.w)  =  -(L(1)Tw,Tw)  Clt  Twll  2 


£  C ( LTw , Tw )  =  C(Tw,w) 

(2.12)  -(TGw,w)  <_  C  (Tw,w)  +  <-  (TGw,  Gw) 

(2.13)  (TGw, Gw)  =  -(G  TGw.w)  <  CllTGwlI  llwll 

<  C  ( LTGw ,  TGw )  1  llwll  =  C(TGw,Gw)^  llwll 
2 

(that  is,  (TGw, Gw)  <_  Cllwll  )  and 

(2.14)  llTfll  =  II  T(L  +  G)  Tfll  £  Cll  Tf  II  . 

Thus,  returning  to  (2.10),  we  can  now  see  that  a  suitable  choice  for  >  0 
leads  to  the  following: 

(2.15)  (Tw,w)  +  llwll2  <  C  (Tw,  w)  +  Cll  Tfll 2  for  0  <  s  <  t  . 

s  —  — 


This  implies  that 


t  t 

/  llwll  2 as  <_  C  J  II Tfll2 as  . 


We  can  now  obtain  (2.5)  from  (2.7),  (2.9)  and  (2.16). 

We  now  turn  to  (2.6).  We  first  see  that  if  0  <  t  <  t,  then 


Ct  llwll  <_  t ( Lw , w ) 

t  t  ...  t 

=  J  (Lw,w)ds  +  /  s(£'  w,w)ds  +2  /  s(Lw,w  ) ds 

0  0  0 


t  t 

<_  C  J  (Lw,w)ds  +  2  /  s(Lw,w  )  ds 
0  0  a 


Since  w  +  Lw  =  f  -  Gw,  we  see  that 


2 

II  w  ||  +  (Lw,w  )  =  (Tf,L*w)  -  (Tf  /  L*w)  -  (Gw#w  )  for  0  <  s  <  i 

s  s  s  s  s  — 

Ihus  integration  gives  us  the  following,  where  >0: 

t  2  t  ,  2 

f  sllw  II  ds  +  2  f  s { Lw , w  )ds  <  —  t(Lw,w)  +  CtllTfll, 

0  s  0  s  -  2  1 


+  C  /  ||  Tfll  ds  +  C  /  (Lw,w)ds  +  c,  f  s  II  Tf  II  ds 

0  1  0  1  o  s  1 


Returning  to  (2.2),  we  find  that 

^-r"l|wl|2  +  (Lw,w)  =  (Tf ,L*w)  <  ClI  Tfll  2  +  i  (Lw,w)  for  0  <  s  <  T 
z  at  —  ll  ~ 

which  gives  us  the  following  estimate: 

t  t  t 

/  (Lw,w)ds  =  /  (Lw,w)ds  <_  C  /  II  Tf ||  ds 
0  0  0  1 

We  can  now  obtain  (2.6),  which  completes  the  proof. 

We  now  study  results  that  hold  for  the  solution  u  of  (1.1). 

Proposition  (2.3):  We  have  the  following  for  0  <  t  <_  T  and  m  >_  0: 

(m+1)  (m)  fm]  (m-J)  (f) 

(2.17)  u  +  Lu  =  -  I  L  u 

9=0  '  '•* 


(2.18) 


Proof :  Since  u  >  C™  ^  (  (0, r  1  , L2  (Cl)  )  ,  we  can  obtain  (2.18)  by  successive  differentia¬ 

tion  of  the  equation  Tu  +  u  =  0.  This  shows  that  u  <  C1"  ( (0 ,  t  ]  ,  D^)  .  Successive 
differentiation  of  (1.1)  then  gives  (2.17)  and  completes  the  proof. 

Proposition  (2.3)  will  be  used  extensively,  in  conjunction  with  elliptic  regularity 
results,  to  translate  information  about  derivatives  in  time  of  u  to  information  about 
spacial  norms  of  u.  For  instance,  an  inductive  argument  shows  that 

™  (? ) 

II  u ( t ) II  f  2m  C  i  ^  u  (till  fol  for  *  1  0  and  0  <  t  <_  1 

1 2m+ 1  j  '=0  (lj 
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We  now  derive  a  series  of  estimates  that  relate  norms  of  the  solution  u  at  one 
time  to  norms  of  equal  or  less  weight  at  an  earlier  time,  giving  up  a  constant  that 
contains  a  pole  in  the  time  increment.  These  estimates  will  prove  important  in  later 
applications . 

Proposition  (2.4)  :  If  0  <  !c  <_  m  and  0  <  s  <  t  £  r,  we  have  that 
(2.19) 


(t  “  s)  ^  li  u  ^  (t)ll  ,  (t  -  s)  II  U  (  t  )li 


2m 


(2.20) 


<  C  l  llu(j)  <s)ll  -  Ctl  u  ( s )  I! 

—  -  2  (m-? ) 

3-0 


(t  -  s)  '  It u (m)  <t)|| 


(t  -  s)  || u  ( t ) | 


2m+  1 


(2.21) 


1C  J  llu<j>  (sHlj  iCllu(s)ll2(m.?)+1 
3=0 


(t  -  s)m+'  l!uim)  (tlllj  1  C|| u  ( s ) II  . 


Proof :  We  will  prove  (2.19)  and  (2.20)  first  and  it  will  suffice  to  prove  these 
results  with  m  >  0  and  i  =  0  or  1,  for  just  the  time  derivatives.  If  •  >  1, 
the  interval  [s,t]  can  be-  split  into  f.  equal  pieces  and  the  results  for  =  1 
used  l  times.  Equations  (2.17)  and  (2.18)  and  elliptic  regularity  can  be  used  to 
obtain  the  results  for  the  spacial  derivatives.  We  note  that  we  are  avoidinq  the  case 
s  =  0,  so  we  have  sufficient  smoothness. 

If  m  >  0  and  0  '■  t  <  r,  (2.17)  shows  that 


1 

2  dt 


u(m)(t)ll"  t  (I.u,m\u'ro')  1C  l  Hub'll.  llu'm'|| 


(m)  (m) 


m- 1 


<3),, 


(m)  „ 


3=0 


m-] 


<  ciiu(m) (tin2  +  c  y  ii u ( j)  (tin 2  . 

3=0 

We  can  now  obtain  (2.19)  for  m  >  0  and  (  =  0  by  an  induction  arqument.  (Note  that 
llu  (till  <  II  u  ( s )  II  for  0  "  s  •  t  -  i). 
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We  also  have  the  following  for  0  <  t  <_  t  : 


(m+1)  „ 2  1  -  (m)  (m) 

l  ( t ) II  +  -  (Lu  ,u  )  t 


m-1 


<  C  l  llu(jll 
j=0 


II  (m+l)..  ,  (tn) ..  2  1  .  (m+1)  2 

-llu  II  +  C||  u  II ,  +  —  II  u  II 
l  14 


<  ~  llu<nl+1,||2  +  C(Lu(m)  ,U(D1) ,  +  cTllu,j,ll?  - 

-  2  j-0  1 


Thus,  we  can  again  use  an  inductive  argument  to  obtain  a  result,  in  this  case  (2.20) 


with  m  >  0  and  1=0. 


(m) 


Now  set  w(t)  E  tu  (t  +  s)  where  m  1 ,  s  >  0  and  0  <_  t  £  t  -  s.  Then 
m-1 


w  +  Lw  =  -t  T 
t  .  „ 

1=0 


L<m-j)u(j>  +  u(m) 


=  -t  "l1  M  -  Y  M  E  f  for 

j=0  j  =  0  >■  1  > 


m-1 

j 


0  <  t  <  T  -  S 


and  w(0)  =  o.  Using  Proposition  (2.2),  we  see  that  if  0  <  £  <  t  -  s,  then 


w(OH  =  C'lu*1"1  (£  +  s) II  <_  C  l  llu®^  (s ) II  +  ~  sup  ;llu(m)  (;  +  s) It 


m-1 


j-0  0<£<£ 

This  gives  (2.19)  with  2  =  1  and  m  ^  1,  after  a  change  of  variables  is  made.  We 
now  derive  (2.20)  the  same  way,  using  the  II  *11  ^ -norm  results  of  Proposition  (2.1). 
To  prove  (2.21),  we  observe  that  (2.17)  and  (2.19)  show  that 


(t  -  S)m+1  Ilu(m)  ( t )  U  2  £  C(t  -  s)m+1  ||L(t)u(m)  (t)  II  <  Cll  u  ( s)ll 


and  that  {t  -  s)  II  u  ^  ( t ) II  <_  C(lu(s)||  ,  for  0  <  s  <  t  <  t.  Thus  we  can  use  an 


interpolation  argument  to  see  that 


(m) 


U . (till  L  1  ca(t)uw(t|,ul"'(t))!  <  Cll  ulm;  (till5  HMt)u'm'(t)ll 


<m>  /mi  i 


/mi  i 


<  Cllu<m)  (till  f  II  u <m)  (till*  <  C  ( t  -  II  u  ( s )  II  . 


This  completes  the  proof. 


We  will  now  study  when  u  is  smooth  at  t  =  0.  The  initial  data  v  must  suM-sfy 
boundary  compatibility  conditions  for  this  to  occur.  We  will  use  certain  (unbounded \ 
"time  dif ferentiation"  operators  A*^  (t) ,  defined  for  m  >  u  and  0  <  t  "  to 

study  these  conditions.  These  operators  will  satisfy  the  equation  u  ^  (tj  -  Afrny  ftiufti 
for  0  <  t  £  t  and  their  form  will  be  motivated  by  (2.17).  Fix  0  t  <  t  ,  let 
(t)  =  I  on  D(A^(t))  =  L^(ft)  and  let 


A""+l>(t)  ;  -  l  H  L<m-3>,t)A!j)(t)  , 

j=0  W 


provided  that  {A^(t)}T_Q  are  given,  be  defined  on  the  domain 

D(A  (t)J  E  {w  :  w  c  D(A^m'(t))  and  A  ^m'(  t)  w  e  D^} 


We  note  that  A^*(t)  =  -L(t),  A^*(t)  =  (L^  -  L^)(t)  and  that  elliptic  regularity 
shows  that  D(A^m*  (t) )  c  H^m  for  »  ^0.  Moreover, 

C  (SI)  5  {f  c  C  (SI)  :  supp  f  c  >;}  c  D(Alm'  (t) )  for  m  >  0 
c  ~ 

We  will  now  use  these  time  derivative  operators  to  characterize  when  the  solution 

u  of  (1.1)  is  smooth  at  time  zero. 

Proposition  (2.5):  (1)  If  v  e  DfA^fO))  and  A^(0)v  e  H*  for  some  m  >  0,  then 

(2.22)  u  €  Cm((0,T] ,H^)  n  C( [0,t) ,H2m+1) 

and  u  e  Cm  ^  ( [0 , T] .D^)  if  m  >  1 . 

(2)  If  v  e  D(A^m*D  (0))  for  some  m  ^  0  then 

m+1  2  m  2m+2 

(2.23)  u  c  C  ([0,t],L  (SI) )  n  C  ([0,t),Dl>  n  C([0,t],H  )  . 

We  note  that  this  Proposition  indicates  when  the  restriction  "s  >  0"  can  be 

removed  in  Proposition  (2.4),  via  the  taking  of  limits,  when  v  satisfies  the  correct 
compatibility  conditions. 

In  the  future,  if  we  say  that  v  e  D(A^m+?*  ( t ) )  for  some  m  ^  0  and  0  _<  t  £  t, 

we  will  mean  that  v  c  DfA^bt))  and  A^m^(t)v  e  H*. 


m  mil  i  * 


Proof :  Our  previous  discussions  indicate  when  equations  of  the  form  (2.2)  have  limits 
in  Dj  as  t  ■+  0;  that  is,  when  wQ  e  and  f  e  C€  ( [0 ,  t!  ,  L2  (;2)  )  for  some  c  >  0. 

Now  we  will  study  the  1*  case.  If  z  c  ,  (2.20)  shows  that  ||U(t,0)z||j  £  C||  z||  ^ 
for  0  <  t  <  t.  Thus,  by  density,  U(t,0)  :  H  -+  is  uniformly  bounded  for 

0  <  t  <  i.  Since  U(t,0)z  z  in  IP  as  t  -  0  if  z  r  D  ,  U(t,0)z  z  in  IP  as 

-  *  L  * 

1  l  •  2 

t  ■+  0  if  z  t  H# .  Thus  if  w  satisfies  (2.2)  where  w  -•  and  f  k  C'([0,t],L  (  )) 

for  some  e  >  0,  w(t)  ■*  w  in  H*  as  t  -*■  0. 

0  * 

Equation  (2.17)  now  suggests  an  induction  argument  which  would  use  the  above 

and  Dr  convergence  results,  which  would  fit  with  our  definition  of  D(A^  (0)) 

L 

and  which  would  complete  the  proof.  However  it  would  be  necessary  to  know  when  the 

2 

right  hand  side  of  (2.17)  is  in  C'  ([0,r],L  (  ') )  for  some  .  •  0,  since  the  induction 

hypothesis  would  only  indicate  that  the  right  hand  side  is  continuous.  The  required 

result  is  discussed  in  Proposition  (2.1),  so  the  induction  argument  can  be  carried  out. 

Thus,  we  can  prove  the  results  on  the  time  derivatives  of  u.  Equations  (2.17) 

and  (2.18)  and  elliptic  regularity  then  complete  the  proof. 

2 

We  now  wish  to  identify  the  !.  (  )-ad joint  of  U,  the  fundamental  solution 
operator.  This  identification  will  prove  useful  Inter  in  some  bootstrapping  nr 'aments 
(as  in  Helfrich  (5)).  Wo  state  it  here  and  sue plv  a  proof. 

Let  L(t)  -  L*  ( r  -  t)  define  another  fainilv  of  differential  operators.  We  note 
that  these  operators  have  all  the  prop*  r*  i«  s  required  of  the  family  fl.(t)  Thus  the 
equations 

U  ( t ,  s )  L(t)P(t,s)  =  o  and  U(s,s)  1  ,  for  0  s  •  t  •-  r 
define  a  fundamental  solution  operator  P  tn  which  parabolic  regularity'  applies.  We 
have  the  following  result. 

Proposition  (2.f>);  b*(t,s)  =  U(  -  s,;  -  t)  for  0  •'  s  t  ^  t  . 

Proof :  Let  U(t,s)  =  II*  f  r  -  s,i  -  t)  for  u  <  s  •'  t.  p  t  and  note  that  —  U(t,s) 
exists  in  B(F.2(  ))  if  t  s.  Thus  if  f  ■  !.“'(  ),  g  «  and  s  •  t,  then 
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(Ufcf,g)  =  ( f ,  U  ( r  -  S,T  -  t)  g)  =  (Uf,L*(t)g)  =  (LUf.g)  ; 

the  last  step  is  valid  since  the  previous  steps  showed  that  Uf  e  D  .  Thus,  U 

satisfies  (2.1) .  Moreover,  if  0<s<t<T  and  f  and  g  e  D  ,  an  estimate  from 

L 

[10]  shows  that  |(U(t,s)f  -  f,g)|  <_  c|t  -  s|  II  f  II  2  II  gll  .  Thus  U  is  strongly  continuous 
on  for  0  '  s  <  t  <  t.  Since  II  Ull  =  Hull  £  1  for  0<_s<^t<^if  U  is  strongly 

continuous  on  L^(:.).  This  completes  the  proof. 

We  conclude  this  Section  with  results  of  a  technical  nature  concerning  the  fine 

structure  of  the  (t) -operators .  Let  0  <_  t  £  t  and  K  ^  0.  Define 

2 

L+(t)  3  L ( t )  +  K  on  Dl  and  let  T+(t)  c  B(L  (fd)  be  the  associated  solution 
operator;  that  is,  L  T  =  I  on  L^(ft).  We  define  L^^(t)  and  tJ^  (t)  as  before, 
for  j  ^  0.  Let  A^(t)  =  I  and,  for  m  0,  define  the  following  (inductively)  on 
D(A(m+1)  (t) )  : 

(2.24)  A<m+1)(t>  -  l  H  =  -  (  l  W  +  KA^>  ,  . 

1=0  (  '  f.=0  ( 

Note  that  A|m) (t) (e  Ktu(t))  =  f^-1  (e  Ktu ( t ) )  if  m  >  0  and  0  <  t  <  t  and 

,,  v  H  2- _  C  2 )  r  fra)  m -<t  Kt  .  -Kt  (1)  (m) 

(2.25)  )  j„j  K  A+  u  =  I  (K  e  )  (e  u)  =  A  u 

if  m  _>  0  and  0  <  t  £  t  .  Letting  v  range  over  D(A*m^  (0))  shows  that  we  can 
take  t  =  0  in  (2.25).  Suitable  translations  of  the  origin  t  =  0  then  show  that 

(2.26)  A ^  (t )  =  2  |  Km  f'A^/  '  (t)  on  D(A  ^  (t) )  ,  for  0  ^  t  £  t 

Now  let  ( t )  i  I,  (t)  ;  L+(t)  and,  for  each  m  >  1,  define  E_^m+^  *  (t)  as 

follows  on  D^,  whenever  each  E^  ^ (t)  r  (e|^ (t) )  *  exists  in  B(L^(Q)),  for 
0  <  j  <  m: 


(2.27)  E  m+U(t)  '  L4(f  +  hit'1  +  l  (-1) 

?=0 

We  have  the  following: 


m+l-f  m  (m-£)  -  (P.+2)  Cm) 

T+  E+  ...  E+  )  . 


Proposition  (2.7):  For  each  m  >_  1,  there  is  a  K  =  K(m)  >_  0  so  that  the-  following 
hold  for  l<^H^_m,  0  <  t  «  i  and  K  ^  K  : 

(2.28)  :  n  ■*  H1  exists  and  is  bounded  for  i  >  0  , 


(2.29) 


-(H) 


•»  H1+2  n  exists  and  is  bounded  for  i  >_  0 


(2.30)  a“  =  (-l)Vm) 


,,(1)  ,  (m) 

E,  on  D(A  ) 


Proof:  Let  0<t<T,K>0  and  f  c  L  (&) .  Since  L  T  f  =  LT  f  +  KT  f, 
— —  —  —  —  +  +  +  + 

(C  +  K)  II T+  fll 2  £  (LT+f,T+f)  +  KllT+fl|2  =  (f,T+f)  <_  II  fll  II T  f  II  . 

Thus  (C  +  K)  ||  T+  fll  Cll  fll.  Also  llT+fl|  <  CllLT+f|l  <  Cll  fll .  Moreover, 

llL+T+ft!  =  II  fll  and  if  l  >  1, 


II L  T<nfll  <  [  (*]  llLU'j)T<j) 

+  +  '3  =  0  W 


fll 


(H) 


so  that  by  induction,  llL+T  II  <_  C  for  2.  >_  0.  Thus 


(?)  f  9  ) 

(C  +  K)  II fll  =  (C  +  K)  11T+L+T^  'fll  <  Cll  fll  for  l  >_  0  . 


We  know  that  is  well  defined,  (2.28)  and  (2.29)  are  satisfied, 

A*1*  =  -E  ^  and  (C  +  K)  llE  ^11  <  C.  We  will  now  assume,  for  some  1  <  1  <  m  -  1, 
+  +  +  —  —  — 

that  E+*^,...,E^^  are  well  defined,  (2.28)  and  (2.29)  are  satisfied, 

A|j)  =  -E^’a^  1(  for  1  <  j  <  (  and  (C  +  K)  II  E~ ( j }  II  <_  C  for  1  1  j  1  1. 

Equation  (2.27)  then  defines  and  the  inductive  assumptions  show  that  (2.28) 


follows  for  E 


U+l) 


Note  that  E^  +  1)  =  L+  +  HL+T.[1)  +  R^  +  1)  =  L+  +  B^+1)  where  (C  +  K)  ilR_['+1)|| 

<  C  and  II L  T(1)ll  <  C.  Thus  we  can  choose  K  >  0  so  that  2llB^  +  ^ll  <  K.  Then 
—  +  +  —  —  +  — 

2  1 
for  any  g  r  L  (.9),  there  is  an  unique  w  e  H  that  satisfies: 


(2.31) 


D(w,se)  +  K (w,t£ }  +  (B+w,v!)  =  (g,^)  for  all  •?  e  H, 
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since  the  form  on  the  left  hand  side  of  (2.31)  is  coercive  over  h|.  But  since 


1 


D(w,^)  =  (g  -  Kw  -  B+w,^)  for  all  >  H4 ,  elliptic  regularity  implies 


that 


w  t  D^.  Moreover, 


(C  +  K/2)  llwll  D(w,w)  +  K(w,w)  +  (B  w,w)  =  (g,w) 


so  (C  +  K)  llwll  <_  ell gll  .  Thus  +  d  ♦L^C':)  is  invertible  and  (C  +  K)  ||E+* 


C.  It  is  now  easily  checked  that  (2.29)  holds  for  E 

U.+ 1) 


(2.30)  . 

On 

,  U> 

+  dT( 

L  A 

+  + 

+ 

1-2 

M-l 

y  (-D 

U+l) 
) ,  we  have  that 


j=0 

=  L  (AU)  +  Y  r 
3=0 


l\  TU’j)A(j  +  1)) 
+  + 


E-U)A(n) 


But  since 


j=o  hJ  +  +  +  >0  hi  +  +  ii0  N  + 


=  T  T  I4]  |j)  lt'H'l'H'a111 
,La  ]  1  +  +  +  + 
L=0  1=1  K  J 


(j 


i) 


i=u  j 

n-i  ( 


-  I 

i=0 

■  i  Til  l+<(t+v 


v‘V  j';1) 

1=0  ^  J  ' 


i=0 
l-l 


-  T  L 
+  + 


)A. 


-  -  I  11 

1=0  1  ’ 


U+l) 


-EU+1)AU> 
+  + 


we  see  that  A+ 

The  proof  can  now  be  completed  by  induction. 


Under  the  conditions  of  this  last  result,  we  see  that  if  m  >_  1 ,  the  operators 
A|m*  can  be  factored  into  operators  that  are  each  a  bounded  perturbation  of  L  and 
which  ,re  in  fact  L  if  we  are  dealing  with  time  independent  coefficients.  Moreover, 
if  we  regard  A^m*  as  an  operator  from  D(A'm')  (to  which  we  give  the  II •  II  ^“norm) 
to  L^(.l)  if  m  ^  0,  it  has  a  bounded  inverse  which  we  will  denote  by  A+^. 
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Ill .  The  Semidiscrete  Approximation 

In  the  last  Section,  we  discussed  some  properties  of  the  parabolic  equation  we 
intend  to  study.  Now  we  are  going  to  examine  a  method  of  constructing  an  approxima¬ 
tion  to  its  solution,  given  a  method  of  approximating  solutions  of  the  associated 
elliptic  problem.  We  begin  by  stating  the  kind  of  properties  we  expect  such  elliptic 
approximation  methods  to  have,  although  we  defer  verification  of  these  properties 
for  several  standard  methods  to  a  later  Section.  We  then  define  the  semidiscrete 

approximation  to  the  solution  of  the  parabolic  problem  and  begin  our  study  of  it. 

2 

We  assume  that  we  are  given  a  finite  dimensional  subspace  S  c  l  (ft)  (that 
depends  on  a  small  parameter  0  <  h  <  1)  and  a  family  of  approximate  elliptic  dif¬ 
ferential  equation  solvers  { T  ( t ) }  that  are  (at  least)  bounded  operators  on 

n  °<tlT 

2  2 

L  (f2).  We  will  give  the  L  (!2)- inner  product.  We  will  further  require  that  the 

following  hold  for  0  <_  t  <  t  : 

(3.1)  (i)  T  <t)  ,  T*(t)  :  L2  (Q)  ■»  Sh  , 

2 

(3.1)  (ii)  (T ,  f,f)  >  0  for  f  •-  L  (Q)  with  strict  inequality  if  0  ^  f  ■  S. 

h  —  n 

Since  L,  ( T,  I  )  1  :  S,  -*  S,  exists,  we  can  ask  that 

h  h  J.  h  h 

h 

(3.1)  (iii)  L*'  1  (t)  -  (^4)  L  (t)  :  S  -*  S  exists  for  1  >  0 

h  'dtJ  h  h  h  - 

Finally,  setting  G,  (t)  3  4  (L,  (t)  -  L*  ( t )  >  ,  we  will  require  that 

n  2  h  h 

(3.1)  (iv)  II G,  (t)^ll  ,  II  ^11  ^  <  C  (L.  (t)^,^),  for  all  0  <  t  <  t  and  <£  *:  S,  , 
h  —  +  n  —  —  h 

where  C+  is  some  (strictly)  positive  constant. 

2  „ 

Given  f  •  L  {-!),  we  will  regard  T^f  as  a  function  m  that  approximates 

the  solution  Tf  of  an  elliptic  differential  equation  problem. 

We  now  make  some  observations  concerning  the  fT^(t)}  and  iL^(t)}  operators. 

2  2 

Let  P  :  L  p)  -»  denote  the  orthogonal  L  P ) -pro jection  onto  and  let 
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Q  :  I 


-  P.  Then  T  =  PT.P,  T*  =  PT*P  and  (3.1)  (iii)  implies  that  the  operator 

n  h  h  h 

T^x^(t)  =  f~*]  (t)  :  L2  <o)  -*  S,  exists  for  any  l  -*  0.  We  let  L.  (t) 

n  ^at'  n  n  —  n 

^  +  Lj*(t))  and  T^(t)  =  (L^(t))  ^P,  for  0  _<  t  <_  x  and  note  that  these 

selfadjoint  operators  are  smooth  in  time  and  satisfy  analogues  of  (3.1).  .Sinc«- 

(L .'fit'?)  =  (L.  'fi,'?)  for  c  S,  ,  we  also  see  that 
h  h  h 


u\\2  <  C((Lj^,s 5)  £  Cll  (L^II  Ml  for  all  y'  Sh  . 

Thus,  II T,  (till  <  C  and  lit  (till  <  C  for  0  <  t  <  t. 
h  —  h  —  —  — 

We  now  use  the  family  { T  ( t ) }  to  define  an  approximation  to  the  solution  u 

of  (1.1).  Choose  a  vh  -  (which  should  be  thought  of  as  an  approximation  to  the 

2  1 

initial  data  function  v  e  L  (ft))  and  let  u^  €  C  ([0,t],S^)  be  the  solution  of 


(3.2)  %  t  +  Lhuh  =  0  £or  0  £.  t  <_  t  and  u^(0)  =  v^  , 

or  equivalently, 

(3.3)  ThUh  t  +  °h  =  °  ^°r  0  i  t  i  T  and  u^tO)  = 

The  function  ( t )  is  our  semidiscrete  approximation  for  u(t).  A  possible  choice 

for  might  be  Pv.  We  will  discuss  other  possibilities  in  a  later  Section. 

We  now  study  some  properties  of  the  solutions  to  equations  like  (3.2).  We  will 

include  estimates  of  the  time  derivatives  of  such  solutions  that  are  independent  of 

the  dimension  of  S,  .  To  enable  us  to  obtain  these  estimates,  we  will  assume, 
h 

throughout  this  Section,  that  the  following  hold  for  0  <_  s,t  <.  x  : 


(s)ll,  II  T  ( t )  l/ ^  ’  ( s )  Plf  <  C  (1)  for  i  >  0  , 

h  h  h  h  —  B 

II G^'  (tMI 2  <  C  «)  (Iyftx,*)  for  t  >  0  and  v'  c  S, 

‘  h  —  B  h  h 

d  i 

where  =  (— )  G^.  Straightforward  calculations  show  that  Condition  then 

holds  for  the  (L^(t)}  family  with,  perhaps,  new  constants. 

We  note  that  (3.1)  and  Condition  imply  that  we  can  use  the  work  of 

Sobolevskii  [10]  to  study  a  fundamental  solution  operator  for  (3.2)  and  Le  assured  o 


dimension  independent  estimates.  Thus,  there  is  a  family  of  operators  U  (t,s) 

h 

V  defined  0  £  s  1  t  <■  t,  that  is  smooth  in  s  and  t  for  all  0  -  s  - 

and  that  satisfies  the  following  for  0<s<(;<t<T: 


(3.4) 


U,  +  L  (t)U  =  0  ,  U.  -  OL  (s)  =  0 

n,t  h  h  h  ,s  hh 

Uh(s,s)  =  I  ,  Uh(t,£)Uh<C,s)  =  Uh(t,s) 


We  will  write  U ^  (t,s)  =  Uh(t,s)  for  0  <  s  £  t  r  T  and  m  >  0.  If 

6  and  Ejj  t  C((O.Tl,Sh),  then  the  following  generalization  of  (3.2): 


(3.5) 


wh,t  +  Lhwh  =  eh  for  0  i  1  £  T  and  wh(0)  “ 


h  ' 


has  an  unique  solution  w  c  C  ((0,t),S  )  given  by 


(3.6) 


Wh(t)  =  Uh(t'0)“h  +  /  Vt,s)V*)ds  ’ 


0 


If  m  0  and  0<_s<t£t,  we  also  have  that 

(3-7)  IIU^1  (t,s)P||  <  C(t  -  s)"m  . 

Versions  of  (2.17)  and  (2.18)  hold  for  u,  and  lead  to  certain  estimates.  For 

n 

instance ,  if  0<t<Tfm>0  and  ^  e  S  ,  we  have  that 

_  _  ft 


.  (m+l) 


.  (m ) 


m-1 


£=0 


(3.8)  (i)  U”’"  w  (t,0)  +  Lh(t)U^‘"'  (t,0)  =  -  l  ™  L^m  °  (t)U^e)  (t,0) 


'  h  ' 


(3.8)  (ii)  II  Lh(t)U^m)  (t.Olv'll  <C  £  II  U  ' 1 ’  ( t  ,OMI 


m+1 


1=0 


(2) 

h 


,  (m) 


We  now  define  operators  (t)  on  s  ,  for  each  0  <  t  <  t  and  m  >  0 , 


that  satisfy  U 


(m)  (m) 


.(0) 


^  U  .  Let  =  I  and,  for  m  0,  let 


A'm+11<t>  -=  -  I  (?j  ^‘(tl^ht) 


(3.9) 


Atlm)(t),'h<t)  +  (Ahm>(t))t 
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{where  the  alternate  characterization  follows  from  the  observation  that  - 

h 

^,t  and  a  soon  to  be  noted  proof  of  the  invertibility  of  U  )  .  We  also  have 
that 


(t)  =  <-I-h  (t)  )m  +  R^m)  (t) 

where  '  (t)  is  a  linear  combination  of  operators  that  are  at  most  an  (n  -  l)-fold 

product  of  ^-ooerators,  includina  at  least  one  time  differentiated  operator. 

If  we  let  L  (t)  t  I  *(t  -  t)  for  0  <  t  \ ,  we  see  that  these  operators 

satisfy  all  the  assumptions  we  made  on  the  family  :L,  (t)},  with  the  same  constants. 

h 

We  can  thus  define  an  associated  fundamental  solution  operator  and  time  differen¬ 
tiation  operators  .  We  also  have  that 

(3.10)  U*(t,s)  =  -  s,T  -  t)  for  0  <  s  <  t  <  i 


(m) 


w  (t) 
h 

let 


We  can  also  use  the 


=  V,  (t,s)w  (s)  for 
n  n 


k,  =  sup  ii  l  mu . 
h  0<-<T  ^ 


energy  techniques  of  Section 
s  <  t  <  :  where  0  *  s  <  t 
Then  since 


II  to  derive  estimates.  Let 

is  fixed  and  w,  (s)  e  S,  and 
h  h 


0  =  T  W  (till  2  +  (L  w  ,w  )  ( t)  •'  (  (  1  -t-  K  ) II  w  !!  2 )  ( t )  for  s  <  t  <  T  , 

2  dt  h  h  h  h  —  2  di  h  h  — 

we  can  conclude  that  U.(t,s)  :  S.  *  S,  exists  for  0  <  s  <  t  <  t.  We  can  also 
h  n  n  —  — 

provide  an  S  -analogue  of  Prooosition  (2.?).  In  fact,  if  vr  e  C1  (( 0,  :]  ,S,  )  solves 
n  p  h 

(3.5)  where  w^(0)  ==  w^  =  0  and  f^  r  C  ([0,r|,S^),  we  have  the  following  for  e  >  0 

(3.11)  II  w  (till  <  C  sup  II  (T.  f,  )  (s)  ||  +  <  sup  II  (T  f  )  ( s )  II  for  0  <  t  <  t  . 
h  —  .  .  h  h  „  hhs  — 

0<s<t  0<s<t 

To  prove  this  result,  we  first  note  that  S  -analogues  of  (2.7)  throuqh  (2.10)  hold 

h 

and  that,  because  of  Condition  B^,  the  following  holds  for  every  >  0  : 


(T,! 1  ’ w.  »w  )  !  =  I  (F,  * 1  ’  T  w  ,T  w  ) 
h  hh  h  h  h  h  h 


(3.12) 


(T.  r.'1):f  w  ,l/1);T w  >  +  C(T  w  ,w  ) 
Ihh  hhh  hh  hhh 


Cc.llw  ||  +  C  (T  w  ,w  ) 

1  h  hhh 
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Then  since  analogues  of  (2.12)  through  (2.16)  hold,  we  can  prove  (3.11).  The  proof  of 


Froposi tion  (2.4)  and  the  invertibil i ty  of  the  operator  than  shows  that  if  rn  ■  0, 


(3.13) 


(t  -  s) '||U1(im)  (t,s)v'H  <  l  llA^j)  (s)v'll 


for  0  <  s  <  t  <  T,  0  <  t  m  and  $  e 

Thus  we  have  (dimension  independent)  estimates  for  the  solution  of  (3.2) 

(and  (3.5)). 

'Cow  let.  K  >  0  and  set  L  (t)  -  L,  (t)  +  K  on  S,  and  T  (t)  =  (L.  (t)  ) 

—  h ,  +  n  h  h ,  +  h ,  + 

on  L **(.'.),  for  0  <  t  <  t.  If  we  define  operators  A^mj  (t)  for  m  _>  0  and 

0  '  t  <  :  using  the  S^-analogue  of  (2.24),  we  see  that  an  analogue  of  (2.25)  holds 

in  P,  .  Also,  the  invertibili ty  of  U,  shows  that  an  analogue  of  (2.26)  holds  in 
h  h 

Finally,  we  can  use  the  techniques  introduced  in  the  proof  of  Proposition  (2.7) 

to  show  that  an  analogue  of  Proi-osi tion  (2.7)  is  valid  in  S.  .  If  m  >  0  and 

h  - 

r.  s  0  is  sufficiently  large,  there  are  invertible  operators  f E,  \  (t)  }m  on  S 

h ,  +  =U  h 

for  0  '  t  <  t ,  given  by  an  appropriate  modification  of  (2.27),  that  satisfy 

(3.14)  A‘m)  (t)  =  (-l)mE‘m)  (t)  ...  E^ht)  for  0  <  t  <•  x  . 

r» ,  +  h,+  h,+  __ 


V  ^  (t)  -  (A^m^(t))  1  and  note  that  IfE.  ^^(t)P|l  <  C  for  0  <  t  <  r. 

h,+  h,+  n ,  +  —  ~  ~ 


We  let 


I 


IV.  Error  Estimates 

We  will  now  study  how  well  the  semidiscrete  approximation  u.  ft)  a:  ;  rex : :  j  ■  • 

the  solution  u(t)  of  (1.1),  given  various  conditions.  W.-  will  continue  *o  '«-n:  .r- 

that.  the  family  (t,  (t)  )  satisfies  (3.1)  but  we  will  no  longer  assume  4  ;ar  ..on  i.t  ; 
h 

necessarily  holds.  We  recall  that  the  somidiscroto  approximation  i,  w.j-: 

defined  in  (3.2)  (or  (3.3)),  where  v,  =  u,  (0)  -  S,  is  thought  of  as  a:,  a;  :  ro>::r-i- 

h  h  h 

tion  for  v.  We  now  let  e(t)  =  u,  (t)  -  u(t)  and  n(t)  "  ( T  -  T,  )  u  ^  ^  (t)  for 

h  h 

0  <  t  <_  t  .  Note  that 

(4.1)  T^e^  +  e  -  p  for  0  <  t  <  t 


We  will  analyze  this  error  equation  in  a  manner  suggested  by  work  in  [3]. 

Our  main  estimate  is  given  by  the  following: 

1  2 

Proposition  (4.1):  Suppose  that  e  and  r  e  C  ({0,t],L  (“.))  satisfy  (4.1)  and 
that  for  any  0  <  6  <  1,  there  is  a  C=C(6)^0  so  that 

(4.2)  j  (Tl[1)  (t)g,g)  |  <  611  g||  2  +  C(6)  (T  (t)g,g)  for  g  t  L^C.)  and  0  <  t  ^  i 
n  —  h  —  — 


Choose 

P  = 

0  or  1  and  if 

p  =  i. 

suppose  that 

Condition  B,  holds.  Then  for 
n 

any 

>  0 

and  0  <  t  <  t# 

we  have 

that 

(4.3) 

tPlle 

(t)l|  <_  c(!  Crfpe)  '0)11  +  C(p 

sup  l|Ci(S)ll  + 

sup  sPllc(s)ll 

o+l. 

+  *  Sur  S*  i'l  .  (  h  )  ■ 

0<s<t  0<s<t 

0<s  t 

where 

o  ( t ) 

t 

•  /  0  (s) ds  for 

0  ^  t 

i  T  • 

0 

We  note  that  (3.12)  shows  that 

Condition  B, 
h 

actual  ly  imj: 

lies  (4.2) . 

Proof : 

Let 

w(s)  -  sm//2e(s) 

where 

m  =  0  if  p  = 

=  0  or  m  = 

3  if  p  =  1.  Then 

(4.4) 

T  w  +  w  = 
h  s 

m/2 

s  a  + 

m  (m-2)/2 

I  S  The 

for  0  <  s 

i  T 

where 

w  (0 ) 

is  e(0)  if  m 

=  0  or 

0  if  ra  =  3 , 

Let  a(t) 

<  C^  ( (0,  x  ]  ,5,  )  be 
h 

defined  by 

(4.5) 

m 

T  a  +  a  =  - 
h  s  2 

(m-2 ) /2 

s  T  e  ;  a  (0)  =  -1 

h 

J’  Pe  (0)  if 
[  0  if 

m  =  0 

m  =  3 
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Thus  if  w  =  a  +  b,  then 


(4.6) 


T  b  +  b  =  sm/20  ;  b(0)  = 
h  s 


Qe ( 0 )  if  m  =  0 
0  i  f  m  =  3 


m/2 


Note  that  s  p  -  b  c  and  that 


(4.7) 


Pb  =  L  (sB1//2p  -  b)  =  L  (sm/2p  -  b)  +  G  (sm^2p  -  b) 
s  h  h  h 


Thus,  if  we  extend  G,  to  an  operator  on  L  (si)  by  the  formula  G  r  G,  P,  we  h 
h  h  n 

that 

(4.8)  T  b  +  b  =  sm^2p  +  sm^2T,  G,  p  -  T  G  b,  for  0  <  s  <  t 

h  s  h  h  h  h 

Let  0  <  t  <_  x  and  <r  >  0.  We  can  use  (4.2)  ,  (4.6)  and  (4.8)  and  the  fact 
that  (T  b,b) (0)  =0  to  show,  as  we  did  in  the  proof  of  Proposition  (2.2),  that 


(4.9)  llb(t)l(2  <_  C(  sup  s1"!!  p  ( s )  II  2  +  sup  sm  2||ps(s)l!2) 

0<s<t  0<s<t 

Thus  if  p  =  0  (so  m  =  0)  ,  the  fact  that  li  a  ( t ) II  <  II  Pe  (0)11  allows  us  to  complete 
the  proof . 

If  p  -  1  (so  m  =  3) ,  we  use  (4.5)  to  see  that 

(4.10)  (T,  a  , a  )  +  ~  -f-  Nall  2  =  ~  s  (m_2) /2  (T.  e ,a  )  for  0  <  s  <  x  . 

h  s  s  2  ds  2  ns  — 

We  now  use  uui  usual  techniques  to  show  that  if  >  0,  then 


m-2 

2 

s 


m-2 
2  * 

(T  e , a  )  =  s  (L.  T  e,T  a  ) 
h  s  h  h  h  s 


m-2 


m-2 


=  s  2  (L  T  e ,T  a  )  -  s  2  (G  T  e,T  a  ) 
hh  hs  hh  hs 


t  (L  T  a  ,T  a  )  +  Cs1"  2(l,T  e,T ,e) 
2  h  h  s  h  s  hhh 


+  cjT  a  II2  +  Csm'2l!G  T  ell2 
2  h  s  h  h 


<  C<  (T  a  ,a  )  +  Csm  2(T,e,e)  for  0  <  s 
—  2  h  s  s  h 


T 
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since  llT^pll  C  (L^T^.T^g1)  =  C(T^,^)  for  >p  e.  and 


IlG.T  tfll  <  C  (L.T  ,  T  ^ )  =  C(T,  •?  ,~f)  for  ~f  <-  S. 
hh  —  hhh  h  h 


Thus,  we  can  show  the  following: 


(4.11) 


11  a { t ) II  <_  C  l  s(c,T  elds 

o  n 


Since  the  error  equation  implies  that 


(T  e  ,T  e)  +  (e,T  e)  =  4  "T-  H  T,  ell  2  +  (e,T  e)  -  (T*^e,f  e) 
hsh  h2dsh  h  hh 


(p  ,T,  e)  for  0  <  s  <  r 
h  — 


( 1 )  2 

and  since  (using  Condition  8,  )  we  have  that  (T,  e,T.e)  <  Cll  T  ell  ,  we  have  the 

h  h  h  —  h 


following: 

(4.12) 


C  ,2,  .  _  2„  ,2. 


f  s  (e  ,T  e)  ds  <  C  f  II T  .ell  ds  +  C  f  s  lip  II  ds  . 

o  h  -  0  h  0 


.(1) 


We  now  note  that  T,e  +  e  *  (Te)  +  e  -  T,  e  =  p  for  0  <  s  <  t,  which 
h  s  hsh  - 


implies  that 


(4.13) 


t  t 

T  e  +  J  e  ds  =  a  +  /  T '  e  ds  +  (T  e)  (0) 


Thus  it  follows  that 


,  .  t  t  t  t  ...  t 

II  rell'1  +  -2-  (/  e,T  /  e)  +  11/  ell  =  (o  +  /  t'  e  +  (T,  e)  (0)  ,  T  e  +  /  e) 

h  dt  0  h  o  0  0  h  h  h  0 


t  t  t 

2  (/  e,  (T  -  T  ) e )  +  (/  e.TV  J  e) 
0  0  n  0 


If  <  >  0,  we  can  use  Condition  B,  to  see  that 

3  h 


(c  +  /  T,(1)e  +  (T.e)  (0)  ,T  e  +  /  e)  <  c  ,  II  Tell  2  +  c  II  /  ell  2 

0  h  h  h  o  "  3  h  3  o 


Cll  oil  +  Cll  (T,  e)  (0)11  2  +  C  /  II  Tell  ds 
h  0  h 
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and  we  can  (again)  use  an  estimate  suggested  by  (2.13)  to  see  that 
t  t 


(/  e,  (T,  -  T  )  e )  =  (/  e,T  G  T  e) 

_  h  n  '  n  h  h 


<  C  ( /  e , T  ( t )  J  e)  +  t  (T  i 
J  r.  h  '  3  h 


W-W’ 


Thus ,  i f 


t  t 

C(f  e /  e)  +  Ct  JIT.  ell  ‘ 
0  h  0  3  h 


is  sufficiently  small,  we  find  that 


(4.14) 


t  ,  t 

/  II T  ell  ds  <  C  /  II oil  ds  +  CtllT.  (0) e  (0) II 
0  h  0  h 


We  now  use  (4.9)  through  (4.14)  to  complete  the  proof. 

Thus  we  see  that  if  we  want  to  estimate  II  u  -  u^ll  under  various  assumptions  on 

v,  it  suffices  to  estimate  quantities  involving  p  =  (T  -  T.)u  and  its  derivatives 

h  t 

We  will  assume  that  the  following  estimates  hold  throughout  the  remainder  of 
this  Section: 


There  is  an  r  >_  2  so  that  if  g  <  H  for  some 
0  <  l  <  r  -  2  and  p  >  0,  then 


(T<P)  (t) 
n 


T  ^  (t) )  gll  <_  (p)  h  +2|l  gll  ^  for  0  <  t  <  t 


If  we  set  p  =  0  in  the  above  inequality,  then  we  are  stating  the  usual  kind 

of  approximation  assumption  for  Galerkin-type  methods.  We  will  show  in  a  later 

Section  that  the  inequalities  are  also  reasonable  for  such  methods  if  p  >  0. 

Our  first  application  of  Proposition  (4.1)  will  be  a  preliminary  result  for 

smooth  and  compatible  data  v.  We  will  need  Condition  B,  ,  described  in  Section  III 

h 

if  we  wish  to  analyze  the  convergence  of  time  derivatives. 

Proposition  (4.2):  Let  m  >  0  and  suppose  that  we  have  one  of  the  following  if  m  = 

(i)  Condition  A  holds  for  the  f T  }  and  {t}  families  and  h  >  0  is 
h  h 

sufficiently  small. 
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j 


(ii)  The  following  estimate  holds  for  some  C*  >  0  : 

(4.15)  i  (L^U  (t)y',vf)  I  <_  C’  I  U^lt)*,*)  I  for  all  0  t  ■_  t  and  r'  . 

or  (iii)  Condition  B,  holds. 

h 

If  m  >  0,  assume  that  Condition  B  holds.  Fix  Q  <_  i.  r  -  2  and  supp  .sv  that 

(a)  f>+2 

v  <-  D(A  (0))  where  a  =  m  +  — — .  Then 

llu'm),t>  -  l^+2^l!+2+2m 

(4. 16) 

+  C  y  llA^hoiv,  -  PA<:i>(0)vil 

j-o  h  h 

Proof :  We  will  use  Proposition  (4.1)  to  obtain  this  result.  We  first  note  that  if 

we  have  Condition  A.  for  the  barred  families,  then 
h 

(T^’g.g)  <  Ch2|lgl|2  +  (T(1)g,g)  <  Ch2llgll2  +  C(Tg,g> 

<_  Ch2||  gll 2  +  C(Thg,g)  for  g  <  L2(Q)  and  0  <_  t  <_  X 
Thus  (4.2)  would  hold  if  h  was  small.  If  (4.15)  holds,  then 

I  (T^’g.g)  |  ic'l  =  C*(Thg,g>  for  g  c  L2(b)  and  0  <  t  : 

Finally,  if  Condition  holds,  (3.12)  implies  (4.2).  Thus  the  hypotheses  of  this 

Proposition  imply  (4.2). 

We  now  observe  that  if  0  <  t  <_  t,  then 

II 0  ( t)ll  <  Ch^^2!!  U(1)  ( t ) II  „  <  Chfc+2|l  u (t )li  „  ,  <  Ch2+2|lv(l.  „  , 

“  P+2  c  +  2 

tllp.(t)ll  <  til  (T(1)  -  T ' 1 1  )  u  II  +  til  (T  -  T.)u  II  <  ChC  +  2|l  vll  ,  . 

t  -  h  t  h  tt  -  i+2 

.  S  +2 

Thus  Proposition  (4.1)  implies  that  lle(t)ll  <_  Ch  11  vll  ^  +  2  +  ^  for  anV 

0<s<_t<^T.  We  can  now  let  s  -►  0  to  obtain  our  result  ror  n  =  0.  Since 
v  f  D^,  we  also  find  that 

11 T  e  11  <  Ch  ^1!  vh  +  dl  v,  -  PvJI  for  0  <  t  <  i 
n  t  —  P.  +2  h  —  — 


I 


We  now  proceed  by  induction.  We  assume  that  m  >  1,  that  v  t-  D(A  (0)) 
i+2 

where  u  =  m  +  — —  and  that 

(4.17)  llThe<P+1)"  <  Chi+2|l  vll  j, +2+2p  +  C  f  II J  (0)  vh  -  PA  ' j  '  (0)  v|l 

for  0  <_  p  <  m  -  1  and  0  <  t  <_  t.  Since  v  e  D(A  ^rn+13  (0) )  ,  we  can  differentiate 
the  error  equation  and  show  for  0  <  t  t ,  that 


(m)  (m)  (m)  (1)  (m)  r  (ml  <m—  j )  (j  +  1) 

T,  e  +  e  =  p  -  mT,  e  )  .  T  e 

h  t  ^  h  >0  (ij  h 


We  now  apply  Proposition  (4.1)  in  the  manner  described  before  and  show  that  the 
following  holds,  for  k  >  0  : 

(4.19)  lle(m)(t)ll  <  Chl+2|lvll  +  Cc  sup  II  To  (m+1 3  ( s )  II 

l+2+2m  0<s<t  h 


+  C  l  llA^3  (0)v  -  PA  3  (0)vll 


<  Chll+2|lvll  sup  II  e  (m)  (s )  II 

i+2+2ra  2 

+  C  J  llA^3)(0)v,  -  PA<3)(0)vll  . 

■  «  h  h 

1=0 

The  estimate  given  by  (4.19)  implies  (4.16)  and  we  can  return  to  (4.18)  to 
complete  the  induction  step.  This  completes  the  proof. 

r 

We  will  now  describe  choices  for  v,  that  will  obtain  0(h  )  convergence  for 

n 

the  error  and  some  of  its  derivatives  if  v  is  sufficiently  smooth  and  compatible. 
The  description  will  be  easiest  if  we  only  wish  to  describe  the  error  and  one  time 
derivative. 

Be foie  we  give  this  first  result,  we  make  an  observation  concerning  the  map 


p,  (t)  -  T,  (t)L(t)  :  d  S.  defined  for  0  <  t  <  t,  often  known  as  the  "elliptic 
1  h  L  h  —  — 

projection.”  By  our  approximation  assumption  A^,  we  know  that 


q  +0 

II  w  -  P,wll  3  II  (T  -  T,  )  Lwll  <  Ch  II  wll  _  for  0  <  t  <  i  , 

1  h  —  P+2  —  — 

P  +2  . 

if  w  H  for  some  0  <  J  <_  r  -  2  and  hence  that 

II  w  -  Pwll  <_  llw  -  P  wll  <_  Ch^  +  ^ll  wll  ^  +  2  for  0  <_  t  <_  \ 

f— ) 

2  ;  ft +  2 

Corollary  (4. 3) :  (1)  Suppose  that  v  e  D (A  (0))(c  H  n  d  )  for  some 

0  <_  P  <_  r  -  2  and  either  (i)  ,  (ii)  or  (iii)  of  Proposition  (4.2)  holds.  Then  if 

P,  +  2 

v,  =  Pv  or  v,  =  P, <0)v  or  v,  *  S,  is  chosen  so  that  II  v,  -  vll  <  Ch  II  vil  . 
h  hi  h  h  h  —  1+2 

we  have  that 

P+9 

(4.20)  llu(t)  -  u^tt)!!  1  Ch  II  ^1  p+2  for  0  —  —  T 

(^ 

(2)  Suppose  that  Condition  holds  and  v  e  D (A  (0))  for  some 

2  2 

0  <  P  <  r  -  2.  Then  if  v  =  P,  (0)v  or  v,  -  T  ( 0 ) L  (0)v  or  v,  e  S.  is  chosen 
—  hi  hh  hh 

P+2 

so  that  IlL^tOlv^  -  L(0)vtl  <_  Ch  II  vHI  j,+4 .  we  have  that 

(4.21)  llu(t)  -  (till  +  Hu(1)  (t)  -  u^1’  (t)ll  1  Ch£  +  2II  vHI  ^  +  4  for  0<  t  M  . 

Proof :  Part  (1)  follows  easily  from  Proposition  (4.2).  As  for  Part  (2),  we  see 
that 

II  v,  -  P  vll  <  CllL.  (0)  (V.  -  P  v)ll  <  CfllL.  (0)v.  -  L  ( 0  >  V II  +  II  (P  -  I)L(O)vll)  . 
hi—  h  hi  —  hh 


The  result  now  easily  follows  from  Proposition  (4.2). 

We  now  wish  to  study  approximation  results  for  higher  time  derivatives.  To 

describe  and  prove  such  results,  we  first  need  to  study  some  properties  of  the 

(O)}^^  and  {A^^(0)}  ^  operators.  We  have  the  following  analogues  of 

Condition  A,. 

h 

Proposition  (4.4):  Let  K  >_  0  and  suppose  that  Condition  holds.  Then,  if 

P 

g  p  H  for  some  0  <_  Z  r  -  2 ,  we  have  that 
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(4.22) (i) 


II  (T^^(t)  -  (t))gll  <_  C  (p) h^+2|l  gll  ^  for  0  <  t  <  t  and  p  >_  0; 

that  is,  we  have  Condition  A,  for  the  {T  (t)}  and  {t,  (t)}  families.  More- 

h  +  h,+ 

over,  if  K  is  sufficiently  large,  we  also  have  that 

(4.  22)  (ii)  II  (E^  ^  (t)  P  -  E+  ^  (t) )  gll  <_  Ch^+^ll  gll  p  for  0  <_  t  ^  t  and  p  1 

Proof:  If  K>_0,  0<_t<^x  and  g  e  H  for  some  0  <  Z  <  r  -  2,  then 


g  =  L+T+g  =  LT+g  +  KT+g.  Thus,  T+g  =  Tg  -  KT  T+g  and  similarly. 


Th,+9  =  V 


(p)  (p) 

KT.T.  g.  Let  p  >  0,  set  E  ^  (T  -  T,  )g  and  suppose  we  know 
h  h,+  —  +  h,+ 


that  II  (t[2>  -  T*^)fll  <_  Cht'ril  fll  if  0  <  j  <  p  -  1  and  f  e  H* .  Then,  sir 


1+2, 


E  =  (T<P)  -  T*P)  )g  +  K  l 
h  1=0 


P-1 


V  (Th,+  -  T+  )g 


-  KT  E 
h 


1=0 


(p-fK  (?) 

-  T  )  T  g 


2 

taking  the  L  (0) -inner  product  of  the  above  with  E  shows  that 

II  Ell  2  <_  C  (K)  hil  +  2|l  gll  e  II  El!  -  K  (T^E,E)  <_  C (K)  h^  +  2|l  gll  p  II  Ell 

We  can  now  prove  (4.22) (i)  by  induction. 

We  will  also  prove  (4.22)  (ii)  by  induction.  We  know  that  (4.22)  (ii)  holds 

with  p  =  1.  We  now  suppose  we  have  (4.22)  (ii)  for  some  p  >_  1  and  all  the 

i, 

intermediate  cases.  Let  0  <  t  <  i  and  g  <  H  for  some  0  <_  l  <_  r  -  2.  If 
K  ?!  0  is  sufficiently  large,  we  can  set  w  E+'P+**g,  w^  E+ *P  *  Pg  and 


P-2 

R  £  l  (-1) 

j=0 


p-j+1  p 

.  j 


L  T(P'3)E-(j+2)  ...  E'(P) 

+  +  +  + 


and  let  be  R‘ s  counterpart  on  S^.  Note  that 


w  =  Tg  -  p  T  ^  w  -  T  Rw  and  w.  =  T,  g  -  p  T*  *  ^  w,  -  T,  R,  w, 

1  +  ++  h  h,+  h,+  h  h,+  hh 


Thus  if  K  is  sufficiently  large,  we  have  that 
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■gift*****- ***■»'>*-* J. r- 


!!  w  -  w  II  <  C  [||  (T  -  T  )  gll  +  ll(T(1)  -  T,(1>)w|| 
n  —  +  n ,  +  +  n ,  + 

+  II  (T  R  -  T  R  P)w|l  +  Ht'1’  (w  -  will  +  II  T.  R,  P  (w  -  w  )  SI  J 

+  h,+  h  h/+  h  h, +  h  n 

<  C(K)h?+2||g||  +  i  || w  -  w  ||  +  ||  (T  R  -  T.  R  P)w||  , 

—  Z  2  n  +  h ,  h 


p-2 


(T  R  -  T.  R  P)  w||  <  C  l  [II  (T* 
+  h,+  n  —  + 

]=0 


<p- j  >  „(p-j),F-(j  +  2) 

h,+  + 


i  Ci)v(rl,p(E-;:)P-E-+,i),A:i>A;^,w„1 


1=3  +  2 


<  C(K)h*+2|lwll  l+2  <  C(K)hll  +  2||gll£  . 

This  gives  (4.22) (ii)  for  p  +  1,  which  gives  the  induction  step. 

We  can  now  complete  our  estimates  for  higher  derivatives  of  the  error  when  v 
is  a  sufficiently  smooth  and  compatible  function. 

Theorem  (4.5):  Suppose  that  m  >_  2 ,  Condition  holds  and  v  e  D(A*a  (0))  where 

o+2 

a  =  m  +  —  and  0  £  Z  £  r  -  2 .  Let  K  >_  0  be  sufficiently  large  and  choose 

w  e  S,  so  that 
h  h 


(4.23) 


"„h  -  A!m)(0,V"  iCh*+2|lvllU2+2m 


for  instance,  let  w  =  PA^m*  (0)v  or  P  ( 0 )  A  ^  (0)  v.  Then  if  v  -  A  ^  (0)w  , 
h  +  1  +  h  h,+  h 

we  have  that 


(4.24)  l  llu(j,(t)  -  u'j)(t)ll  <  ChU2|lvll 

j-0 


£+2+2m 


for  0  <  t  <  t 


Proof:  Proposition  (4.2)  and  (2.26)  and  its  S  -analogue  show  that  it  suffices  to 

-  h 


know  that 


,Ah!X!“>wh  ■  pAij)\<m)AIra)vii  ich*+2nvii 


£+2+2m  ' 


for  each  0  £  j  <  m  -  1.  But  (4.22)  shows  that  the  following  holds  for  0  £  j  <_  m  -  1 


and  g  £  H  : 
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.A^X.^rg  -  A'i,A;(m,qll  <  ch?+2|lgll  .  . 

Since  II  a/  ,  A,  ^Pll  <  C  for  0  <  j  <  m  -  1  and 

n,+  h,+  —  —  — 

II  (P  -  I)A(:i)v||  <  Ch*+2||v|l  ^  _  for  0  <  j  <  m  -  1  , 

+  '  £+2+2:  -  - 

we  have  our  result. 

We  now  will  discuss  convergence  results  for  the  error  and  its  derivatives  when 

2 

v  is  no  more  than  a  function  in  L  (ft).  We  will  limit  our  choice  of  initial  data 

for  the  semidiscrete  approximation  to  vL  =  Pv  for  these  nonsmooth  data  results. 

h 

Theorem  (4.5):  Suppose  that  Condition  B,  holds.  Then  if  v  =  Pv  and  m  >  0.  wc* 
-  h  h  — 

have  that 

(4.25)  y  llu<j)(t)  -  u‘j)(t)ll  <  Chrt_r/2_m!l  vil  for  0  <  t  '  i  . 

j-o  h 

Proof :  We  begin  by  noting  that  the  estimates  of  Sections  II  and  III  and  density 

00 

arguments  allow  us  to  assume  that  v  e  C  (ft) .  These  estimates  also  show  that  we  can 

c 

2 

assume  that  h  <_  t. 

We  now  use  Proposition  (4.1)  with  p  =  1.  Since 


I  a  (till  =  11/  (T  -  t)  (s ) u  (s)dsll 


<  II  (T.  -  T)u(t)ll  +  II  (T  -  T)  u  (0)1!  +  J  ||  (t/1)  -  T(1))ullds 
-  h  h  0  h 

2 

^  Ch  llvll  for  0  ^  t  <  t  , 

2  2 

and  similarly,  tllp(t)ll  +  t  <_  Ch  II  vil  for  0  t  <_  i ,  we  see  we  have  shown 


lu(t)  -  u,  ( t ) II  =  II  (U  t,0)  -  U.  (t,0)P)vll  <  Ch2t  *11  vll 
h  h  — 


for  0  *  t  <  t.  Thus  we  have  shown  (4.25)  for  m  =  0  and  r  =  2.  To  obtain  (4.25)  with 
m  _>  1 ,  we  will  need  to  use  this  result  and  our  previous  smooth  data  results  in  a 
bootstrapping  argument.  We  will  also  use  a  special  representation  for  smooth  and 


-33- 


compatible  functions.  (This  representation  was  used  in  a  slightly  different  context 


,(P) 


in  ID).  Suppose  that  p  1  and  w  t-  D(A  *  (t))  for  some  0  <  t  m.  Then,  if 

K  >  0  is  sufficiently  large, 


(4.27) 


-d) 


■  (!)„,  „(!> 


w  =  A  ,r'PA'r'w  +  (E  '  '  -  E  P )  E  w 

h,+  +  +  h,+  + 


-  I  AT(i-1,P(E-<i)  -  E-(^P,A(j)w 
,L.  h,+  +  h,+  + 

3=2 


Let  0^p£m,  0<h  <_  t  £  T,  t  =  t/3,  t  =  2t/3  and  v  e  •  Then 


lu(P'  (t)  -  u<P’  (till  =  II  (U!P>  (t  ,0)  -  u‘P>  (t  ,0)  P)  vll 
n  n 


£  II  (U<p)  (t,t2>  -  U<P)  (t,t2>P)U(t2,0)vll 


+  llu'p)  (t.t.)Plt  [IUU(t,,t.)  -  U  <t,,t  )P)U(t  ,0)v)|| 
h  2  21  n  2  1  1 


IIU(t_,t  )  (U  (t,  ,0)  -  U.(t,  ,0)P)vll 
2  1  1  n  1 


+  !lu(t,,t, )  -  U.  (t,,t,)P||  llult,  ,0)  -  U.  (t.,0)Pll  II  vll]  . 
21  h  2  1  1  hi 


-p 

<_  term^  +  Ct  (term2  +  term^  +  term^) 


By  our  smooth  data  results  applied  with  a  suitable  translation  of  time  zero, 


we  find  that  if  K  is  large,  then 


term  <  II  U (p) u (t  >  -  U *p)  aT  <P>  PA jP)  u  (t  )ll 
1—  2  h  h,  +  +  2 


y  llu'P  a”  i"*'P(E.  -  E,  '^'P)A'J'u(t  , 

L  hh,+  +  h ,  +  +  2 


j-i 


ICh  llu(t2)llr+2p 


tc  [  t'p~1+jhr||u(t,)ll 


3  =  1 


2  r-2+2 j 


<  Chrt'r/2'Pllvll  . 


r  - r/2 

We  observed  in  (4.20)  that  term^  <  Ch  t  II  vll  .  Also, 


1 

3 


term  <  II U  (U  -  U,P>I|  ||  vll  =  II  (U*  -  U;P)U*II  II  vll 
i  n  h 


=  II  (0  -  0,  p > 011  II  vll  <  Chrt"r/2||vll 
n  — 


where  we  have  used  our  adjoint  identifications  and  smooth  data  results  applied  to 
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the  time-reversed  operators.  Finally,  term^  £  C(h  t  )  II  vll  by  (4.26). 


By  iterating  the  above  argument  u  times  with  p  =  0  we  find  that 


?  “1  r /7  7  -1  7 

IU(t,0)  -  U.  (t  ,0)  Pll  <  C  ( (h  t  )  '  +  (h  t  )  )  , 

h  ~ 


.U+l 


for  any  0  •<  t  <_  t  .  We  now  choose  y  >_  1  so  that  2  >  r.  The  proof  is  then  oasil 

completed  by  red'  -ng  the  argument  once  for  each  1  <_  p  £  m. 

We  conclude  this  Section  by  noting  the  following  result  on  forcing  terms.  If 

we  assume  we  are  given  a  suitably  smooth  function  u(x,t)  on  Cl  *  [0,1]  that  satisfies 


u^  +  L(t)u  =  f(t)  on  2  *[0,t]  and  u(0)  =  v  on  SI  , 


where  f  is  suitably  smooth  on  Q*(0,t],  we  can  define  a  semidiscrete  approximation 
by  the  following: 


u,  +  L  u  =  Pf  for  0  <  t  <  t  and  u,  (0)  =  v, 
h ,  t  h  h  —  h  h 


where  v  is  chosen  from  S,  .  (We  note  that  u,  (t)  always  exists).  If  we  set 
n  n  h 


e  :  u^  ■  u,  we  find  that 


T,  e  +  e  =  (T  -  T.  )  (u  -  f)  =  -  (T  -  T)  Lu  for  0  <  t  <  t  . 
h  t  h  t  h  — 


This  equation  can  be  easily  analyzed  by  Proposition  (4-1).  For  instance,  if  we  set 


v,  =  Pv  and  assume  Condition  B,  ,  then 
h  n 


I  u ( t )  -  u  ( t ) II  <  C(u)h 
h  ~ 


for  some  constant  C(u)  depending  on  the  solution. 


V.  Examples 


We  will  summarize  here  some  well  known  results  ...  v  r  :!  i. 

projection  methods  and  we  will  give  the  additional  require.!  irti.tf-  t'-r 
methods  that  will  allow  us  to  apply  the  theory  of  the  preceding  . 

by  sketching  some  of  the  common  features  of  the  methods. 

Each  method  will  use  a  finite  dimensional  subspace  of  functions  m  i;‘ 
will  be  associated  with  parameters  0  <  h  <  1  and  an  r  >  2  in  the  following 

(5.1)  min 

*  c  sh 


»ll  +  hllw 


Ch£+2llwl! 


;  +  2 


f  +  2  1 

for  all  wen  r  H  ,  where  0  <  £  <  r  -  2  and  where  II  >11  is  a  (perhaps  h- 

—  —  I 

dependent!  seminorm  that  will  be  related  to  the  IMI  ^  norm  but  which  may  contain 

other  terms  dealing  with  boundary  condition  considerations. 

For  each  of  our  coercive  differential  operators  L(t),  there  will  be  an 

associated  positive  form  □.  (t) (*,*)  that  is  related  to  the  Dirichlet  form  of  the 

h 

operator  and  which  will  be  used  to  define  the  associated  T,  (t)  operator.  Given 

h 

2 

f  ■  I.  ,  the  function  w^  =  T^ltlf  t  will  be  defined  by  the  following 
equations : 

D.ttXw  ,v-)  =  (f,v-)  for  V"  e  S.  . 
h  h  h 

Tin-  form  D  (t>  will  be  symmetric  positive  definite  if  L  =  L*  =  L. 

The  following  relations  will  hold,  with  certain  constants: 


(5.1  ) 

cplvll  2  ■_  M2  ± 

C.D  (t)  (y',v 
2  h 

)  for  all  ^  c  S.  , 

h 

(5.4) 

!Dhl,(t,(9i'92>! 

1  C3  OOllg^l 

ill  g^H  j  for  Py  >_  0  , 

where  g. 

and 

g2  t  +  Dl  and 

D'e,(t,(... 

)  denotes  the  form  obtained 

by 

tak  ing 

successive  time  derivatives  of  the 

D,  -form.  We  will  also  have 
h 

the 

following 

for 

j  >  0  : 
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|  (L(j)w,g)  -  Dt*  ■'  *  (w,  7) 


(g.L*(j,w)  -  r^1 


where  w  ..  H  n  for  some  ■"> 


£  +  and  z  •  D  ir 

h  L  L 


arbitrary.  (Thus  we  will  "almost"  bo  able  to  integrate  by  parts),  finally,  wo  will 
have  that 

d  1  ? 

(5.6)  G  (t  )~f  a  P  (  V  a  .  (t )  D.  y  +  —  (  ,  3.  a  .  (t ) )y  )  for  c  <•  S  , 

h  01  1  2  1  Ox  h 

2  2 

where  P  :  L  (2)  -*■  is  the  usual  L  (.'.) -orthogonal  projection. 

Vte  will  now  list  vhat  the  II  *11  -  semi  norm  and  the  D  -form  is  for  each  method. 

I  h 

We  will  then  go  on  to  prove  the  necessary  stability  and  accuracy  results.  (There 
are  references  given  in  [2]  and  [3]  for  the  various  methods.  We  will  not  repeat 
them  here ) . 

(1)  Gale rk in's  Method  for  the  Neumann  Problem: 

In  this  case  D  (t )  =  D(  t),  II')!  =  ]J  •))  and  C.  can  be  taken  to 

h  I  1  * 


(2)  Gale rk in's  Method  for  the  Dirichlet  Problem: 

Now  we  must  restrict  to  be  in  Thus  boundary  conditions 


are  required  of  functions  in  S  .  Again  D  (t)  =  D  (t ) , 


I  =  II  •  II  and 


(3)  A  Method  of  Nitsche; 

This  is  a  technique  for  the  Dirichlet  problem  that  does  not  require 
that  lie  in  H^.  The  norm  and  forms  for  the  method  are  as  follows: 

D.  <t )  (- ,  - )  £  D(t)  (•,•)  -  (•)>  -  <4~  (•),•)  +  .-h~h  •  ,  •  )  , 

h  jn  3n_ 

11  ’ " i  -  »■»!  +  cih"1|l'llo,ifi  +  ciih""-)n"o,:.:  +  IMIi,-.  1  - 

a 

where  r —  denotes  the  conormal  derivative  associated  with  L  (t )  ,  {•) 

3n  n 
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denotes  the  normal  derivative,  C  and  C  are  certain  constants  and 

3  >_  0  is  sufficiently  large.  The  following  inverse  assumption  is  required 

on  S.  as  well: 
h 

11  +  lkl1.3fl  -CIIl(h'i|U"ill  +  h_1H  v"  q  <  all  *  -  Sh  . 

Wa  note  that  we  can  take  =  0. 

Another  Method  of  Mitsche : 

This  is  another  method  that  handles  the  Dirichlet  problem  without 
requiring  that  lie  in  H^.  The  method  is  the  same  as  the  previous 

one  but  3  may  be  taken  to  be  zero.  However  the  following  "almost  zero 
boundary  conditions"  restriction  must  be  put  on  : 

IMI  „  <  C  h-IMI  ,  for  all  •f  e  S 

U/iu  —  iv  l  h 

where  C  is  sufficiently  small. 

A  Lagrange  Multiplier  Method  of  Babuska: 

This  is  yet  another  approach  for  approximating  the  solution  of  the 

Dirichlet  problem  without  imposing  boundary  conditions  on  functions  in 

Sh-  The  space  is  constructed  in  a  special  way  so  as  to  agree  with  a 

Lagrange  Multiplier  formulation  that  would  be  used  in  practice.  Wfe  will 

not  detail  this  construction  here,  but  note  that  the  key  to  the  boundary 

f +2  1 

conditions  is  given  by  the  following  estimates.  If  w  e  H  n  Hq  where 
0  ^  r  -  2 ,  g  ►:  S  +  and  z  a  is  arbitrary,  then 


.n.D.w,g>|  =  j(  }'  a.  .n.D.w,  g  -  z>| 


Ch  II  wfl  ,  +2 II 9  -  zl!  for  1  <_  i  <_  d  , 

due  to  certain  boundary  aoproximation  properties.  Vfe  use  D,  (t )  =  D(t) 

h 

and  l!*ll  =  11*11.  for  this  method  and  (5.5)  follows  from  (5.7). 
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We  will  now  prove  estimates  that  will  hold  for  any  approximation  method  that  is 


defined  by  (5.2)  where  (5.1)  and  (5.3)  through  (5.6)  hold.  Vte  first  see  that 
(5.8)  |  l  fjjl  Dh!j_l)  (t)  ((T^1’  (t)  -  T(l)  (t ) )  f  ,^2  ) !  ^  CCjl*2  -  2||  Ihi+1|l£llJ  , 


(5.9)  (L^j)  (t)*lW>2)  =  (t)  felf^2)  «  (^i'Lh<j)  (t)^2>  ' 

£ 

for  j>0,0<t<T,v;1  •'f-,  c  S.  /  z  £  Dr  and  f  E  H  ,  for  some  0  <_  2  <_  r  -  2. 
—  ~  1  2  h  L 

The  following  is  our  main  approximation  result. 


Theorem 

(5.1): 

For  m  >  0,  0  <  2  ; 

(5.10) 

||  (T(m)  (t)  -  T!m)  (1 
n 

(5.11) 

||  (T<m)  (t)  -  Tto)  C 
h 

Proof : 

W5  will 

first  prove  (5.10) 

(5.11). 

Say  0 

<  l  <  r  -  2  and 

for  0  < 

j  <_  m. 

bet  e  S.  b< 

h  h 

when  we 

replace 

w  by  w«\  for 

Set  v(j)  E  T(j)f  and  w 'j  ’  E  T'j)f 

h  h 


(m)  (m)|,2  /v  (m)  (m)  v  (m)  (m) 

cll*h  -wh  "iiV#h  ~wh  ’  wh  -wh  > 

<  a(fc6,)  -  w(m),  fc<m)  -  w'm))  +  CCjlw'm)  -  w‘m)llhU1|lfll. 
—  hh  h  h  *h  hi  x. 


1 1  (,] 

j=0 


<  ciife(m)  -  wu(m)n  (chSL+1iifii „  +  y  ii w' 

—  h  hi  2  h 


<j)  -  w(j)ll  ) 


This  can  be  used  to  show  (5.10)  for  all  m  >_  0. 

,  _  *  (m)  (m) 

Now  we  use  a  duality  argument.  Let  z  e  D  satisfy  L  z  =  w  -  w  and 

L  n 

let  i  r  S  be  the  function  that  attains  the  minimum  in  (5.1)  when  we  replace  w 
h  h 

by  z .  Then 
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(m)  (m).,  2  (m)  (m) 

II  w  -  w,  II  =  (w  -  w^  ,L*z) 
h  n 


"  im)  (w(3>  -  „<3>,L*(m-i>z)  -  "j1  [m]  (w<3>  . 

=n  UJ  h  (3 J  h 


1  1 


N  r>-3>  (w(3)  _  „(j) 


j  =  0 


( 3  J  h 


‘  WhJ,'Z  "  V  +  CC*  L  llw<j)  ■WhJ)  !h  •’ 

3  =  0 


L 

j=0 


m-1 


ml  (m-3)  (3)  3  .  .  ~  .,  3  (]),,  ,, 

.  D,  (w  -  w,  ,z,  )  +  C  /  II  w  -  w,  II  II  zl 

3 1  h  h  h  .  h 

j=o 


m-1 


,„,2+2..  ..  _  r  11  (3)  (j  )„  (m)  (m),, 

(Ch  II  fll  „  +  C  l  II  w  -  w  II  )H  w  —  w  II 

j=0  h  h 


This  can  be  used  to  complete  the  proof. 

we  also  have  the  following: 

Proposition  (5.2):  If  i.  >_  0  and  0  ^  t  M  vhen 

(5.12) 


IlG^’ttMl2  <  C  (L  (t)y  r'r  )  =  C  (L,  (t)'r  ,y  )  for  all  'r  t  S. 
h  —  n  n  h 


Proof:  Vfe  see  that  if  *  *  S.  , 


II  (t)*l|2  <_  dlv'll2  idly'll2  <_  CtL^tJy'.y'J 


we  will  now  study  the  remainder  of  Condition  B,  ,  using  inverse  properties 

h 


-1„ 


(as  was  also  done  in  [7]). 

Proposition  (5.3):  Suppose  that 

(5.13) 

Then  for  1.  >  0  and  0<^s,  t  <_  T ,  we  have  that 

(5.14) 


lly'll  <  Ch  lly'll  for  all  #  £  S. 
I  —  n 


LU)  (t)Tu(s)ll  ,  II T.  (t)LU)  (s)Pll  <  C(t)  . 
n  h  h  h  — 


2 

Proof:  If  g  c  L  (3),  ^  c  S,  and  l  >  0,  then 

-  h  — 

|  (L^-)  (t)Th<s)g,*>)  [  =  |D^'  ’  (t)  ((T  (s)  -T(s))g,y')|  +  |  D^0  (t)  (T(s)g.v') 

<  Cll  (T,  (s )  -  T(s))gll  JMI_  +  |  (L(f'!  (t)T(s)g,y')  1  +  CJill T (s)gll _ll y'll 
—  n  II  z  i 

<_  Cllgll  lly'll  . 
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This  suffices  to  prove  the  first  part  of  (5.14)  and  the  rest  of  the  proof  follows 
by  considering  the  adjoint  problem. 

-2 

W6  note  that  if  we  have  (5.13),  llL^(t)Pll  <_  Ch  for  0  <  t  < 

Thus  we  can  now  apply  the  analysis  of  Sections  II  through  IV  to  many  Galerkin 
type  projection  methods. 


-41- 


VI .  Maximum  Norm  Estimates 

2 

We  will  now  examine  how  we  can  use  our  L  -based  estimates  in  conjunction 
with  maximum  norm  estimates  for  the  associated  elliptic  problem  to  prove  similar 
maximum  norm  estimates  for  the  parabolic  problem.  We  will  use  techniques  similar  to 
those  introduced  in  13]. 

We  first  study  global  maximum  norm  estimates.  We  assume,  as  usual,  that  we  have 
a  family  of  operators  -T^(t)  ’•  that  satisfies  the  properties  listed  in  (3.1)  and 

the  approximation  assumptions  .  We  will  also  assume  that  we  are  working  in 

5  - 

d  =  1,  2  or  3  sr-ace  dimensions  and  that  S,  c  C"  {  ?.)  for  some  constant  0  <  6  <  1; 

h  — 

that  is,  the  functions  in  5,  are  Hdlder  continuous  with  exponent  5  on  ft. 

h 

Finally,  we  will  assume  the  following  inverse  property  on  : 


(6.1) 


llv'l 


CU)h'd/2_£llv'  II, 


for  0  <  e  <  4  and  e  S. 


C 


To  obtain  maximum  norm  estimates  for  the  parabolic  problem,  we  need  to  know  some 
corresponding  estimates  for  the  elliptic  problem.  We  will  assume  the  following: 


There  is  a  function 
we  have  that 


Y  (h) 
P 


so  that  if  T(t)f  e 


(6.2)  { 


II  (T  -  T)  (t)  f  II  <  Y  (h)ll  T  (t)  f  II 
h  0 , 00  "*  p  P 


for  some  0  <_  t  £  t, 


|  where  p  =  2  or  r 

We  will  also  assume  that  y^  (h)  <_  Ch ^  ^  for  some  n  <  “•  Work  in  f4]  or  ]9]  done 

2  r 

under  various  conditions  suggests  that  we  could  take  y  (h)  =  Ch  ,  Yr(h)  -  Ch  if 

r  >  2  and  Y_(h)  =  Y  (h)  =  cjloq  h(h2  if  r  =  2. 

2  r 

We  note  that  these  conditions  imply  the  following: 


IIT.v'IL  -  II  (T.  -  T)*|l  +  llTcll 

h  0 ,  —  h  0 ,  *  0,w* 

•-  v  (h)IITrfll  +  Cll-Nll  -  y  (h)IUII  ,  +  Cllc 

2  2.”  ‘  C  1  (0.) 

-3/2-, 

'  (C>  (h)h  <  Oll.'ll  ClL-  II 


whe  re  0 


<  min 


1 

2 


n,  ;-)  and 


S 


h 
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If  we  define  u.(t)  t  S,  via  (3.2)  with  some  choice  of  v,  e  S.  as  usual,  then 
h  h  h  h 

(4.1)  shows  that  if  0  <  t  <_  r , 

(6.3)  llu(t)  -  u.  (t)|l  <  II  (T  -  T,  ) u  II  +  IIT  Pe(1>ll0 

h  0,“  —  h  t  0,°°  h  0,“ 

<  y  (h) (1  u(t)ll  +  Clle(1)  (till  . 

—  r  r 

We  can  now  use  the  results  of  Section  IV  to  further  analyze  (6.3).  For  instance, 

if  we  have  Condition  B,  and  we  set  v,  =  Pv,  we  find  that 
h  h 

(6.4)  llu(t)  -  uh(t)ll  <C(t  )(y  (h)  +  hr)llvll  . 

for  0  <  tg  <_  t  i .  Other  results  can  also  be  formulated  for  sufficiently  smooth 
and  compatible  v. 

Similar  estimates  can  be  done  in  the  interior  of  SI  if  the  appropriate  estimates 
are  known  for  the  { ( t ) }  family. 
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